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The effects of gauge interactions in graphene have been analyzed up to now in terms of effective 

models of Dirac fermions. However, in several cases lattice effects play an important role and 

^^ , need to be taken consistently into account. In this paper we introduce and analyze a lattice gauge 

Ch ■ theory model for graphene, which describes tight binding electrons hopping on the honeycomb 

lattice and interacting with a three-dimensional quantum (7(1) gauge field. We perform an exact 

Renormalization Group analysis, which leads to a renormalized expansion that is finite at all orders. 

The flow of the effective parameters is controlled thanks to Ward Identities and a careful analysis of 

^^ ' the discrete lattice symmetry properties of the model. We show that the Fermi velocity increases up 

Cn I to the speed of light and Lorentz invariance spontaneously emerges in the infrared. The interaction 

produces critical exponents in the response functions; this removes the degeneracy present in the 

non interacting case and allow us to identify the dominant excitations. Finally we add mass terms 

D ' to the Hamiltonian and derive by a variational argument the correspondent gap equations, which 

; , , have an anomalous non-BCS form, due to the non trivial effects of the interaction. 
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1. INTRODUCTION 

Graphene, a two dimensional crystal of carbon atoms that has been recently experimentally realized [45, 46], has 
highly unusual electronic properties. The lattice in graphene has a honeycomb shape and the corresponding energy 
bands intersect at two points, close to which the effective dispersion relation is approximately conical and similar to 
a "relativistic" one. The low energy excitations of the half-filled system consist of hole-particle pairs created close to 
the tips of these two cones. These quasi-particles behave like two-dimensional (2D) massless Dirac fermions: for this 
reason the infrared (IR) properties of the system can be understood in terms of a model of 2D Dirac particles in the 
continuum [51, 56]. As a result, a number of concepts and features of high energy physics have a correspondence and 
can be observed at much lower energies in this crystal. 

The description in terms of Dirac fermions is quite accurate in the free case, and it is also helpful in the presence of 
many-body interactions (see, e.g., [7] for a review), as it allows to translate and adopt a number of powerful methods 
from the realm of quantum field theory (QFT) to that of graphene. However, taking the effective description too 
seriously has some drawbacks, since a model of interacting 2D Dirac fermions in the continuum has spurious ultraviolet 
divergences due to the linear bands. In order to make the continuum theory well-defined, ad hoc regularizations must 
be introduced to cure the short distance singularities, which arc obviously absent in the tight binding model, where 
the honeycomb lattice acts naturally as an ultraviolet (UV) cut-off. It is unfortunate that the computation within 
the Dirac model of certain physical observables, such as the conductivity, is sensitive to the specific choice of the 
regularization scheme, a fact that makes the comparison with experiments difficult or inaccurate, see, e.g., [35]. 

These ambiguities make an approximation-free analysis of the effects of the lattice and of the non-linear bands 
in graphene highly desirable. In this paper, we study a lattice gauge model for graphene that describes electrons 
hopping on the honeycomb lattice and weakly interacting with a three-dimensional (3D) quantum U{1) gauge field. 
Our model has two independent parameters: the bare Fermi velocity v and the electric charge e (we use units such 
that the reduced Planck constant h and the speed of light c are equal to 1). The analysis is performed by using exact 
Renormalization Group (RG) methods, which allow us to express the physical observables in terms of renormalized 
expansions in the electric charge with finite coefficients at all orders, uniformly in the volume and in the temperature 
(more precisely, we show that the coefficient multiplying e^" in the renormalized expansion grows at most as n\). Our 
main physical predictions are the following. 

1. Thanks to the validity of exact lattice Ward Identities, the gauge field remains massless and the IR behavior 
of the system is characterized by a line of fixed points (i.e., the effective charge has vanishing beta function). 
Correspondingly, the physical parameters are strongly renormalized by the interaction. In particular, the Fermi 
velocity increases up to the speed of light and Lorentz invariance spontaneously emerges in the infrared. This is 
proved by fully taking into account the discrete lattice symmetries, which are used to exclude the presence of 
dangerous extra marginal or relevant terms in the RG flow. 



2. The wave function renormalization diverges at the Fermi points with an anomalous exponent. This last properties 
strongly resembles one of the crucial features of one-dimensional Luttinger liquids; in this sense, the model 
considered in this paper is one of the very few established examples of Luttinger liquid behavior in two dimensions 
(it has been suggested that also the Hubbard model on the square lattice close to half filling is a Luttinger liquid, 
but this is still an unproven fact). 

3. The response functions have an anomalous behavior expressed in terms of non trivial scaling exponents. The 
difference between the interacting and non-interacting exponents is small at small coupling. In particular, the 
response functions associated to fermionic bilinears, which decay as r~^ at large distances in the non-interacting 
case, remain integrable even in the presence of weak interactions; therefore, magnetic, phonon or superconducting 
susceptibilities are finite and no evidence for quantum instabilities is found, in agreement with the fact that the 
fixed point is close to the trivial one at weak coupling. 

4. On the other hand, the interaction removes the degeneracy in the decay exponents of the response functions: 
some exponent increase and some other decrease and this depresses or enhances the effects of local perturbations 
associated to specific fermionic bilinears. This gives us a criterium to extrapolate the qualitative behavior of 
the system at larger values of the electric charge and decide what are the favorite quantum instabilities at 
intermediate to strong coupling. An explicit computation of the anomalous exponents shows that the dominant 
excitations correspond to: (i) Kekule distortions, associated to a dimerized Peierls' pattern, (ii) charge- density 
waves associated to an excess/deficit of the electron density on the two sublattices of the honeycomb lattice, 
(iii) Neel antiferromagnetism, (iv) the Haldane circulating currents [27]. In all these cases, the logarithmic 
singularity at zero transferred momentum in the first derivative of the response function in momentum space 
is changed into a power law singularity with an anomalous exponent. The singularity in the other responses is 
either absent or weaker than in the four cases mentioned above. 

5. If we add a symmetry breaking field coupled to a Kekule distortion the induced energy gap in the spectrum 
is dramatically amplified by the interaction: the ratio between the energy gap and the amplitude Aq of the 
external field diverges as Ag — > with an anomalous power law. 

6. The effect of the electronic repulsion on the Peierls-Kekule instability, usually neglected, is evaluated by de- 
riving an exact non-BCS gap equation, from which evidence is found that the gauge interaction facilitates the 
spontaneous distortion of the lattice and the gap generation. Similar conclusions can be drawn for the mass 
terms corresponding to the other dominant excitations. 

The exact Wilsonian RG methods we use are based on ideas borrowed from constructive QFT and have been 
introduced in the context of interacting Fermi systems in [3] (and then extended in [47, 54]) at the beginning of 
the 90s, and since then successfully applied to various problems in solid state physics, see, e.g., [39] for an updated 
review. In particular, they appear to be very well suited to analyze the properties of graphene without any Dirac 
approximation, large-N approximations or unphysical regularization schemes. Such methods have been first applied in 
[15] to the Hubbard model on the honeycomb lattice, where we proved in a full non-perturbative fashion the analyticity 
of the ground state of the half-filled system; we used the same methods to rigorously establish the universality of 
the optical conductivity of graphene with weak short range interactions [20, 21], an issue that still needs to be fully 
understood in the case of electromagnetic interactions. Technically, the analysis in this paper extends the one in 
[18], the main differences being that here we fully exploit the lattice symmetries and, besides computing the reduced 
density matrices, we construct the response functions, we compute the critical exponents and discuss the effects of 
external symmetry breaking fields. Most of the results in this paper were announced in [19]. 

Let us add a comment on the range of applicability of our theory. Our analysis is based on resummations of pertur- 
bation theory in a, where a = e^ /{Airehv) is the effective fine structure constant of graphene, which unfortunately is 
not small: e.g., in suspended graphene, a ~ 2, which makes graphene an intrinsically strongly coupled problem, apriori 
not accessible to approaches based on power series expansions. However, one needs to take into account that the 
effective Fermi velocity is considerably increased by the interactions (recent experiments [12] can observe a factor three 
amplification of the velocity close to the Fermi points and even a larger enhancement is expected at lower energies), an 
effect that goes in the direction of decreasing the effective fine structure constant. Therefore, our theory is valid close 
to the IR fixed point, with effective parameters e and v {v close to the speed of light) that should be thought as being 
obtained by the (non-pcrturbative) integration of the "first few IR scales" , possibly by using numerical methods, like 
those of [9, 10]. Let us also remark that since we predict the emergence of anomalous critical exponents, our final 
results can be easily extrapolated to intermediate coupling, which would not be the case if the apparent logarithmic 
divergences emerging in perturbation theory were not correctly resummed at all orders. 

The rest of the paper is organized as follows. In Section 2, we define the model, we present in detail our main results 
and compare them with existing literature. In Sections 3 to 9 we discuss the proof of our results: in Section 3 we 



discuss the functional integral representation of our model and derive the relevant Ward Identities (Wis) ; in Section 
4 we describe our RG scheme to compute the functional integral; in Section 5 we show how to use Wis to prove the 
vanishing of the photon mass, the vanishing of the beta function for the effective charge, and how to control the flow 
of the renormalization parameters (Fermi velocity, wave function renormalization, vertex functions) ; in Sections 6 and 
7 we compute the response functions associated to several fermionic bilinears; in Sections 8 and 9 we explain how 
to modify the general RG scheme to include the effects of an external field coupled to local order parameters and 
how to derive the non-BCS equation for the gap. A number of more technical aspects of the proof are deferred to 
the Appendices: in Appendix A we derive the functional integral representation, we derive the Wis and we explicitly 
show the equivalence between the Feynman and Goulomb gauges; in Appendices B and G we analyze the symmetry 
properties of the theory and use them to identify the symmetry structure of the relevant and marginal kernels; in 
Appendix D we perform the lowest order computations of the beta function and of the critical exponents; finally, for 
completeness, in Appendix E we check at lowest order the cancellation of the photon mass and of the charge beta 
function, which follows from the general gauge invariance of the model. 



2. A LATTICE GAUGE THEORY FOR GRAPHENE 

In this section we introduce the model, define the main quantities of interest and present our main results. A 
comparison with existing literature is also presented. 



A. The model 

We let A — {ni^i + n-il-i : rii ~ 0, . . . ,L ~ 1} be a periodic triangular lattice of period L, with basis vectors 
li — ^(3, a/S), h = 5(3, — \/3). We denote by Aa — A and A^ = A + ^^ the A- and B- sublattices of the honeycomb 
lattice, with 5i the nearest neighbors vectors defined as: 

li = (l,0), 52 = ^(-l,V3), 53^^{-l,-V3). (2.1) 

We introduce creation and annihilation fermionic operators for electron sitting at the sites of the A- and B- sublattices 
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FIG. 1. The honeycomb lattice of graphene. 
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where Bl = {k = riiGi/L + n2G2/L : < n^ < i}, with Gi_2 = ^(Ij iv^), is the first Brillouin zone; note that in 
the thermodynamic limit L'^ Y^keBL '^ '^'^^ fe ^^' ^'^^ B = {k = ^iGi + C2G2 : C* G [0, 1)} and \B\ = 87rV(3V3). 



The operators a- ^, 6- ^ satisfy the canonical anticommutation rules, and are periodic over A; their Fourier transforms 
arc normalized in such a way that, if /c, fc' G Bl, 

{a+ ,a+ } = {aZ ,aZ } = {bi ,b+ } = {bz ,bz } = 0, {aZ , at } = fbz ,6+ } = L^Sr i:,S^ a> . 

(2.3) 
Definition Eq.(2.2) implies that a- , bi are periodic over the reciprocal lattice A*. 

We also introduce a quantized photon field living in the 3D continuum. Let Sl,L' = Sl x [^\, "y) ^'^ ^ subset of 
R^ with periodic boundary conditions and Sl = {x = LS,\lx + ^^2^2 : 6 •= [0, 1)}; let also Vlx' = T^l x |t^ be the 
corresponding dual momentum space, with V^ ~ \j> ~ niGi/L + niGijL : Ui G Z} (the honeycomb lattice can be 
thought as being contained in the section 5l x at z = of Sl)- For all p — {p,P3) G Vl^l' we introduce bosonic 
creation and annihilation operators c^^, with helicity index r = 1, 2; they satisfy the commutation relation 

[c+^,C+ ,,,] = [Cp^r,Cp, ,,,] = , [Cp,r,C^'y] = L'\SL\Sp,p,6r,r' ■ (2.4) 

Let A{x) ~ {A{x) , Arj{x)) be the quantized vector potential on Sl defined as: 



^(-) = I7^ E Ey^^p-(s,.^"^''^ + sV"^)' (2-5) 

where £p^r G ^^ are polarization vectors satisfying the conditions 

which reflect the choice of the Coulomb gauge. Note that, in the thermodynamic limit |5l|^^ l^nev "^ (^tt)"^ L3 dp. 
Moreover, the function x(|p|) acts as an UV cutoff function: more exactly x(t) is a smooth compact support function 
equal to 1 for < t < | and to for t > 1. The photon UV cutoff is chosen to be on the same scale as the inverse 
lattice spacing. We expect that this cutoff should be removable, but this is not our main concern here. 

The interacting electron-photon system we are interested in is described at half filling and in the Coulomb gauge 
by the following grandcanonical Hamiltonian: 

Ha = H^ + h{ + Va, (2.7) 

where the first term is the (gauge-invariant) hopping term, the second represents the field energy, and the third is the 
Coulomb interaction, namely: 



H'. = -t E E 



f ief^ld.sAix+.sSjfi)-S.i + .- , -ief^d.sAix+.s5jfi)-5j,+ - \ (q R) 
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where t > is the hopping strength, e is the electric charge, and ng is equal to X](T=ti ^s a'^s a "-"^ ^^ 'l2a=i-i ^s a^s a 
depending on whether x G A/i, A^, respectively; moreover, ip{x) is a regularized periodic version of the 3D Coulomb 
potential: 

where we remind the reader that p ~ (pjp^). Note that the electron-photon interaction is induced both by the complex 
hopping rate texp{ie/g dsA{x + sSj) ■ Sj} in Hj^ and by the static Coulomb interaction Va; the combination of the 
two describes the retarded electromagnetic interaction mediated by 3D photons. If the electric charge e is 0, then the 
Hamiltonian decouples into a sum of two quadratic terms: 

^aL.o = -^ E E (4^+.-.. + ^-+,;, .«.-..) +17^^ E \p\4..rC;..r=--Hl + Hi, (2.10) 



which can be expHcitly diagonahzed; the corresponding correlation functions can be computed exactly, via the Wick 
rule, in terms of the electron and photon propagators, which read as follows. Let ^J o- — i^i a^ ^t.+s a) ^^^ ^x cr = 

1-^'^ I be row and column spinors, with ^i — (0, (5i), x ~ {xo,x), xq G [0,/3) an imaginary time, /? > the 
inverse temperature, and a^ „ — e^'^^^a^ q-Haxq 5± = ^HAXofj± ^ ^-H^xa ^j^^^ imaginary time evolved of the 

creation/annihilation operators. Then the /ree electron propagator is [15] 

^o^-^(x) := (T{*,,,*+ J) I =-^ J2 ^"^'^^o(k) , 5o(k) := ^—— ( % """J^^)' 

' ^"^ " ^ kee^ kl + v^\n{k)\^ \-vn{k) tko 

(2.11) 

where (•) „ ^ denotes the average with respect to e~^^A(*)^ x is the fermionic time ordering, Bfs^L '■= %"(^ + ^) ^ ^l, 

w = |i is the bare Fermi velocity and fi(fc) = 2^ gife((5j-iSi) ^^le complex dispersion relation. The function fl{k) 



is vanishing if and only if fc = jyf,, where pyf, — f 4|^, ±^j= I arc the two Fermi points, close to which rt{k' + j?^) = 
ik'i ± fc2 + 0(|/c'p). Therefore, setting p"^ = (0,p^), k' = k — p"^, cj = ±, the propagator in momentum space reads 

where Z = 1 is the 6are wave function renormalization. Eq.(2.12) has the form of the propagator for massless Dirac 
fermions in 2 + 1 dimensions. 

Similarly, defining, for x ~ (xq, x), A^ = e^A^" A(x, 0) e^^A^" , the in-plane free photon propagator is, for i,j G {1,2}, 



«..--(x):=(T[(A.MAo),]),^,L^„ = ^ E E e---|^ %-p^ , (2.13) 



where V/s^l -^ ^^ x Vl- In the limit L' ^ cx) 



»5'(x) = T^^ T. e--"4'"(p) ■ -f (p) - / l^i^K - ul?;;.) . (2-") 



where the apex (C) reminds the choice of the Coulomb gauge. Note that the IR singularity of the in-plane photon 
propagator in momentum space is ~ |p|~^, rather than the usual ^ |p|~^ of standard quantum electrodynamics 
(QED). Therefore, interacting graphene at low energies is similar to a gas of massless 2D Dirac particles interacting 
via a modified |p|~^ photon propagator. 

B. Response functions 

Our goal is to understand the behavior of the system in the presence of a non-zero electron-photon coupling. We will 
be mainly concerned with the computation of the interacting correlations, in particular of the interacting electronic 
propagator and response functions. The latter are particularly relevant from a physical point of view, since we can 
read from their long distance behavior (or, equivalently, from their singularities in momentum space) the tendency of 
the system to develop quantum instabilities associated to several putative local order parameters, in the same spirit 
as [52]. For illustrative purposes, we restrict our attention to the response functions associated to the following bond 
fermionic bilinears: 

C^, - E. (e^^ ^« '' ''^-^''^ <aK+s, ,. + c-c.) (lattice distortion) 

C^f "^ = E. (<.«x,. - bt+s, ,aK+s, ,a) (staggered density) 

C^J ^J2a'^ (aiaflx^a - bt+s^ .ab^+s, .a) (staggered magnetization) 

C^ = Ea f<a«x,. + fex+6,,.fex+6,,a) (bond density) (2-15) 

Cxj = E<T i^e"'^" '^''*''^"+°'''ax,^^x+5,,a + c-c.) (bond current) 

^x,j — l^cr y(^ ■'"x,cr"x+mj,<j *e "x+Sj,a"x+Sj+mj,<j ^ '-■'-■J 

(Haldane circulating currents) 



where in the last hne mi = ^2 — <^3, ni2 = <^3 — ^i and 1113 

corresponding response functions are defined as: 



^1 — S2 indicate next to nearest neighbor vectors. The 



?(a) 



i?i;^(x-y) =. (C,;C;,,) - lim lim {Q,;Q, 



fS^ooL- 



•■J' I3,L 



(2.16) 



where the semicolon in (•;•) indicates truncated expectation: {A; B) := (AB) — {A){B). From the long distance 
behavior of Ry^' (x) we can read the possible emergence of long range order. For instance, if the lattice distortion 

response function behaved as Rlj (x) ~ Ccos(p^(x — Si + 5j)) for some C ^ asymptotically as |x| — > 00, this would 
signal the spontaneous emergence of a Peierls' instability in the form of the dimerized Kekule pattern of Fig. 2a, which 
is one of the possible distortion patterns of graphene [13, 33]. In fact, C^- is the local order parameter coupled in the 
Hamiltonian to the hopping strength, whose mean value heuristically represents the intensity of the hopping from a 
given site to its nearest neighbor. Having C ^ means that far away bonds are strongly correlated and that the joint 
distribution of their hopping rates is non uniform, but rather oscillating with a cosine dependence. With reference 
to Fig. 2a, this oscillation can be heuristically understood by associating a factor proportional to 1 (to — i) to the 
double (single) bonds and by averaging over the three equiprobable configurations obtained by rotating Fig. 2a by 
0,^,^. In momentum space, this would correspond to the appearance of a delta singularity (i.e., of a "condensate") 

in Rlj '{p) at p = p^. Similarly, a condensate in the p = mode of R-- (p), ^L (p) or R-- '{p) would signal 

the spontaneous emergence of a staggered density pattern {charge density wave), of a staggered magnetization pattern 
{Neel order) as in Fig. 2b and Fig. 2c, respectively, or of the specific pattern of circulating currents discussed in [27]. 
Even in the absence of condensation, that is of delta-like singularities in i?j-^'(p), the possible loss of regularity in 

Ry^ (p) due to the interaction can be interpreted as a tendency of the system to develop quasi- long range order in 
the corresponding channel. 




FIG. 2. a) The Kekule distortion; the double and single bonds correspond to higher and smaller hopping rates, respectively, b) 
The charge density wave instability; big dots correspond to a charge excess, while small dots correspond to a charge deficit, c) 
The antiferromagnetic instability; the arrows represent the spins of the electrons, and they have to be understood as lying on 
the axis orthogonal to the honeycomb lattice. 



C. Critical behavior and anomalous exponents 

We are now ready to state our main results. As shown in the following, after systematic resummations of per- 
turbation theory, we are able to express the observables of our theory as renormalized series in the electric charge, 
with finite (and explicitly bounded) coefficients at all orders; the coefficient of e^" grows at most as n!, a behavior 
compatible with Borel summability of the theory at weak enough coupling strength. 

In particular, in the limit j3,L — ?> cx), the interacting two-point function in the Feynman gauge (see Section 3) is 
given by: 



S{x) = (T(V'-,V'o+ J) = 
1 



dkn 



dk 



-zpx 



^(k) 



5(k'+p^) = 



Z(k' 



> (27r) Je \B\ 

iko v{k'){-ik[ + ujk!^) 

v{]i'){ik[ + ujk2) iko 



(1 + S(k')) 



(2.17) 



where ^(k') and w(k') are the interacting wave function renormalization and Fermi velocity, while -B(k) is a sub- 
dominant term, vanishing at the Fermi points k' = 0. The interacting propagator close to the Fermi points has a 
structure very reminiscent of the free propagator, Eq.(2.12). However, Z(k') and w(k') are strongly renormalized by 
the interaction: 



^=T^2+0{e'), ^= +0(e4), (2.19) 



where 



z(k')^ik'r\ i-«(k')=^(i-«)ikr, (2.18) 

where 

p2 r, 2 

liTT-^ OTT^ 

are anomalous critical exponents, well defined at all orders in renormalized perturbation theory (the O(e^) remainders 
in Eq.(2.19) are written below in terms of a series in e^", with coefficients growing at most as n!). Eqs.(2.17)-(2.19) 
are very reminiscent of the IR behavior of the 2-point function of a Luttingcr liquid [28, 41, 52], and are consistent 
with results obtained in a model of interacting Dirac fcrmions in the continuum [22] . 

On the contrary, the interacting 2-point function for the photon has the same IR singularity ^ |p|^^ as the free 
case; this means that the gauge field remains massless (no screening). 

Regarding the response functions associated to fermionic bilinears, we find that in general the interaction changes 
their decay exponents at large distances as compared to the non-interacting case, where all the responses decay as 
^ r^^ at large distances. The presence of the interaction makes these exponents non trivial functions of e (anomalous 
dimensions). In particular we prove that 

<"""^ W - l^^-w^^pi^ + .r"' W > (2-21) 

^5r w - l^^-^TR^ + -i""' w ' (2.22) 

^^""^ - ^ + 0{e'); e(^^'^) = |J + 0(e^); ^^^^^ - ^ + 0{e'); e"^^^ + 0{e^) (2.24) 

and A^ = A^{v,e) arc constants that arc equal to 1 at the free Dirac point, i.e., A^\ =1 (in particular, for 

V close to 1, A^ = 1 + 0(1 — v) + 0(e^)). Moreover, the correction terms r^^ (x) are subdominant contributions, 

decaying at infinity faster than |x|^^+'> ; they include both the effects coming from the irrelevant terms in a RG 
sense and the effects proportional to 1 — w(k') (see Eq.(2.18)) coming from the Lorcntz symmetry breaking terms. 
From Eqs.(2.20)-(2.22), we see that the decay of the interacting responses in the K, CDW, AF, H channels is slower 
than the corresponding non-interacting functions; i.e., the responses to K, CDW, AF, H are strongly enhanced by the 
interaction. On the contrary, all other responses decay at infinity faster than |x|~'*+('=°"''')^ , i.e., if a = D,J, (and 
similar bounds are valid for other observables like the Cooper pairs, see Section 7 below) 

l4"H-)l ^ ^^ ' (2-25) 

for some constant C > 0. These results can be naturally extrapolated to larger values of the electric charge, in which 
case they suggest that the lattice distortion, staggered density, staggered magnetic order and the "Haldane circulating 
currents" are the dominant quantum instabilities at intermediate to strong coupling strength. 

As we noticed, quantum instabilities are also signaled by divergences in the Fourier transform of the response 
functions. Even in the presence of non trivial exponent, the power law decay remains integrable at weak coupling: 
therefore, no divergence is found in the Fourier transform of the response function. On the other hand, 

dpRlf\p' + p±) ^ IpT^*"' and dpR^\p) ^ |pr«'"' , with a = CDW,AF,H , (2.26) 

so that the singularity in the first derivative of the response functions in momentum space, which appeared as a 
discontinuity or at most as a logarithmic divergence in the non interacting case, is enhanced and turned into a power 
law singularity by the interaction for these four responses. The singularity at different momenta or for other response 
functions is either weaker or absent. Also in this respect, the conclusion is that the system shows a tendency towards 
Kekule, charge density wave, Neel ordering or to the formation of the Haldane gap. 



D. Mass renormalization and gap equation 

The enhancement of the response functions suggests that the effects of small external staggered fields coupled to 
the K, CDW, AF, H local order parameters are dramatically enhanced by the interactions. This is in fact the case. 
Let us, for instance, add an external staggered field coupled to the lattice distortion local order parameter, with the 
same cosine dependence as the long distance decay of R\j ^(x), see Eq.(2.20). Physically, this can be interpreted as 
a fixed distortion of the lattice into a Kekule pattern as in Fig. 2a. In fact, if we allow distortions of the honeycomb 
lattice, the hopping becomes a function of the bond length i^j that, for small deformations, can be approximated 
by the linear function tg^j = t + g{(-3.j ~ £) ='■ t + (t>s,j, where i is the equilibrium length of the bonds and (j)3j 
plays the role of a classical phonon field. If a Kekule distortion of amplitude proportional to Ao/g is present, then 
0K,j =32 cos{pp{x — Sj + Sj„)) for some jo G {1, 2, 3} and the Hamiltonian becomes 

^a"=^a-^ E ^ e*P--(^"-*-^+^~^o)C.^,- . (2.27) 

j = l,2,3 xeA 

If the electron-photon coupling e is equal to 0, then the fermionic 2-point function has a mass proportional to Aq. 
If we switch on the interaction, the mass (i.e., the decay constant describing the exponential decay of the 2-point 
function at large distances) becomes 

A = Aj/(i+''-) , r^^^^ + ..., (2.28) 

that is, the Kekule mass is strongly amplified by the interaction (note that the ratio between the interacting and bare 
masses diverges with an anomalous exponent as Aq — > 0). This phenomenon is very reminiscent of the spontaneous 
mass generation phenomenon in QFT, the main difference being that while in a truly relativistic theory in the 
continuum the flow of the mass can be studied also in the UV region and the bare mass can be let to zero with the 
UV cutoff (still keeping the same dressed mass at fixed IR scale), see [25], here the lattice acts as a fixed UV cut-off, 
so that the mass is amplified but not spontaneously generated. Similar arguments can be repeated for the other 
dominant excitations. 

Let us finally discuss a possible mechanism for the spontaneous generation of a Kekule instability in our model. 
Rather then fixing the distortion pattern (j)gj once and for all, we can let 4> = {'/'SjIseA ' be a classical field to be 
fixed self-consistently, in such a way that the total energy in the Born-Oppenheimer approximation is minimal, i.e.. 



= argnnn-^ 



{^o(^) + i E 4.}' (2-29) 



xeA 

i=l,2,3 



4K) 



where Eq{4>) is the ground state energy of H'^ = i?A — X) Seh <Px jCg , • We find that the Kekule distortion pattern 

- J = l,2,3 ' '"' 

4'? = 00 + ^ Ao COS {p + [x - 5, + 4 )) (2.30) 

is a stationary point of the total energy, provided that i^o — c^g^ j k + • • • for a suitable constant cq and that Aq 
satisfies the following non-BCS gap equation: 

« 7 fc2+i;2(k')|0(fc'+p^)|2 + |A(k')|2 

A<|k'|<l 

where A = t^^^^^'^'^^ and, for A < |k'| < 1, Z(k') ~ |k'|-'', w(k') - 1 - (1 - w)|k'|'' and A(k') - Ao |k'|-''^. Our 
gap equation has the same qualitative properties of the simpler equation: 

1 = 5'/ dp—^^ -^ (2.32) 

from which it is apparent that at small e, the equation admits a non trivial solution only for g larger than a critical 
coupling gc] remarkably, gc ~ ^/v, with v the free Fermi velocity, even though the effective Fermi velocity tends to the 
speed of light. Therefore, at weak coupling, the prediction for gc is qualitatively the same as in the free case [33]; this 
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can be easily checked by noting that the Fermi velocity w(k') is sensibly different from v only for momentum scales 
exponentially small in v/e^. See Section 9 for more comments about this point. 

Even more interestingly, the value of ^c decreases as e increases; i.e., interactions facilitate the formation of a Kekule 
pattern. Eqs.(2.31)-(2.32) can be naturally extrapolated to intermediate coupling: if in such a regime rj/s^ — rj — 

j^ + ■ ■ ■ exceeds 1, then the integrand in the r.h.s. of the gap equation diverges as A —t' 0, a fact that guarantees the 
existence of a non-trivial solution for arbitrarily small g. In other words, the larger the electron-photon interaction, 
the easier is to form a Kekule patterned state; it is even possible that at intermediate coupling gc = 0, which would 
imply a spontaneous generation of the Peierls'-Kekule instability. Note the non BCS-like form of the gap, similar to 
the one appearing in certain Luttinger superconductors [40] . A similar analysis can be repeated for the gap generated 
by the staggered density, the magnetization or the Haldane mass. 

E. A comparison with existing literature 

Before we enter the technical part of our work, let us conclude this expository section by a comparison of our model 
and our results with existing literature. We do not pretend to give a full account of the rapidly expanding literature 
on the effects of interactions in graphene; several excellent reviews already exists, like [7, 38], which we refer to for 
extensive bibliography. Here we focus on the difference between the approaches and results based on the effective 
models of Dirac gas in the continuum, which is the most popular and widely studied model of graphene, and ours, 
which is based on a tight binding lattice model. 

Short range interactions. Graphene with electron-electron screened interactions has been studied in terms of an 
effective model of 2D massless Dirac fermions interacting with a local quartic potential: in the weak coupling regime 
this interaction is irrelevant in the Renormalization Group sense [24], while at strong coupling analyses based on 
large A'^ expansions found some evidence for quantum critical points [29, 30, 49, 50, 53]. This model requires an UV 
regularization and some observables, like the conductivity, appear to be sensitive to the specific UV regularization 
scheme used: different results are found [31, 32] depending on whether momentum or dimensional regularizations is 
chosen. 

A more realistic model for graphene with short range interactions is a tight binding model that keeps into full 
account lattice effects, such as the half-filled Hubbard model on the honeycomb lattice. Formally, the Hubbard model 
reduces to the continuum Dirac gas in the limit as the lattice spacing goes to zero; in this sense the latter can be 
thought as a scaling limit approximation of the former. One important advantage of the lattice model as compared 
to the continuum one is that within the former no ambiguities arise in the computation of the conductivity. In 
particular, all the interaction correction to the conductivity exactly cancel out in the optical limit [20], in agreement 
with experimental results, [43], and in disagreement with the Dirac model with momentum regularization. 

Gauge-invariant electromagnetic interactions. In the early paper [22] (written much before the actual realization 
of graphene) an effective model for interacting of graphene was proposed, in which massless Dirac fermions in the 
2D continuum are coupled to a quantum 3D photon field, with the fermionic propagation speed much smaller than 
the speed of light. The main result of [22], based on second order perturbation theory (and, therefore, valid in the 
weak coupling regime) , is that at low energies the Fermi velocity tends to the speed of light and the wave function 
renormalization diverges with an anomalous exponent. No computation of the response function exponents was 
performed. 

In [18], we revisited the model proposed in [22] and, rather than using dimensional regularization as in [22], we used 
an UV momentum cut-off; this is a much more natural choice, since the Dirac continuum model should be though of 
as an effective model emerging in a Wilsonian RG after the integration of the high energy degrees of freedom. Using 
this model, we extended the results in [22] at all orders, but still we computed neither the response functions nor the 
gap equation. On the other hand, the momentum cutoffs, necessary to avoid spurious UV divergences, break gauge 
invariance and this imposed the introduction of counterterms in order to keep the photon mass vanishing and in order 
to have one (rather than three) effective charge [18]. 

In the more realistic model considered in this paper the electrons live on a honeycomb lattice and interact with a 
quantum photon field living in the 3D continuum. The fact that gauge invariance is not broken has the effect that 
no unphysical counterterms need to be introduced. The critical exponents and the gap equation have been computed 
here for the first time. Moreover, lattice gauge invariance prevents the generation of several potentially dangerous 
marginal and relevant terms. As in the case of Hubbard interactions, the conductivity computed in the continuum 
model show an unphysical dependence from the conductivity [35], while considering the model and the formalism 
introduced in this model will resolve such ambiguities. 

Static Coulomb interactions. The most popular model used to describe graphene with unscreened electromagnetic 
interactions is a Dirac gas with static density-density Coulomb interactions; this is an apparently sensible approxi- 
mation, since the (bare) propagation speed of quasi-particles in graphene is much smaller than the speed of light, so 
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that retardation effects should be negligible at least in a wide range of energy scales. However, in the weak coupling 
regime, a second order RG analysis predicts an unbounded growth of the Fermi velocity in the IR and the vanishing of 
the effective charge at the Fermi points [23]. Therefore, at low energy scales the Fermi velocity becomes comparable 
with the speed of light and the model with static interactions loses its significance, a fact that can be seen as a dra- 
matic manifestation of its "uncompletness" , see [55] . This is also consistent with the fact that the theory with static 
Coulomb interactions docs not appear to be rcnormalizable at all orders, see [19, p. 1425]. In conclusion, retardation 
effects are important to understand the nature of the IR fixed point of the theory. 

The effective Dirac model with Coulomb interactions has been also extensively analyzed in the strong coupling 
regime. It has been argued that, at large enough coupling, an excitonic gap spontaneously opens [25, 36, 37], by 
a mechanism similar to mass generation in QED2+1 [42]: in these works, the gap equation is derived by a self- 
consistence argument; the corresponding solution is shown to be momentum-dependent and vanishing in the limit of 
high momenta. However, these findings rely on several approximations, in particular: in [25, 36] the vertex, wave 
function and velocity renormalizations are neglected; in [37] the renormalization of the velocity is taken into account, 
but the corresponding flow is IR-unbounded and UV-cutoff dependent. In our work, the gap equation is obtained 
by using an exact energy optimization problem and fully takes into account the renormalization of all relevant and 
marginal operators; moreover, since we do not neglect effects of the honeycomb lattice, our results are free from the 
ambiguities related to the presence of the spurious UV divergences typical of the Dirac approximation. 

In [26, 30] a systematic classifications of the possible interaction and mass terms allowed by symmetry is performed 
in the Dirac model with Coulomb interactions; in the present paper a similar analysis is carried out, with the difference 
that the lattice discrete symmetries rather than the continuum symmetries are taken into account. 

RG analyses based on large N expansions [11, 29, 30, 53] and Quantum Montecarlo analyses [9, 10] in the presence of 
Coulomb interactions have identified critical exponents for the response functions and found evidence for the presence 
of excitonic phase transitions (like CDW and Kekule instabilities). Again, these results are in qualitative agreement 
with our finding that the effective Kekule mass (or the CDW, AF, H mass) grows at low momenta with an anomalous 
power law, although the model and the method are quite different (expansion in the charge and retarded interactions 
versus 1/N expansions and instantaneous interactions). A strong coupling expansion for a lattice gauge theory for 
graphene has been also performed in [1, 2] and evidence for spontaneous Kekule mass generation was found. It is 
worth stressing that the lattice used for the Quantum Montecarlo analysis in [9, 10] is a square lattice, rather than 
the original honeycomb lattice; it would be interesting to repeat a similar analysis for the more realistic honeycomb 
lattice gauge theory introduced in the present paper. 

Let us finally note that the Peierls-Kekule instability was first discussed in the non- interacting case in [33], as a 
key ingredient for the emergence of electron fractionalization without the breaking of time-reversal symmetry (on this 
issue, see also [8, 34]). Our gap equation generalizes the one of [33] to the interacting case. 



3. FUNCTIONAL INTEGRAL REPRESENTATION AND WARD IDENTITIES 

The correlation functions introduced in the previous section can be conveniently expressed in terms of a functional 
integral. We introduce the generating functional of correlations in the ^-gauge with infrared cutoff on the photon 
propagator as 

gW«'''*(*,J,A) ^ /"p(-^^)p€./i*(-^^)gV(*,A+.7)+e(*,A+J,<I.) + (A,*) ^ (-3 -^^ 

where: 

1. ^ = {'I'xo-} ^^'^ two-components Grassmann fields, the components being denoted by 'i'^ ^ „, p = 1,2, the first 
corresponding to the a^-fields, the second to the 6*-fields. Moreover, A — {A^^^}, with fi — 0, 1,2, are real 
fields. The convention on the Fourier transform of the fields that we use is the following: 



.T.± \ ^ ±ikx,f, A 



^ E e"'^%.P- (3-2) 



2. If n^ = 6f,fi, ^ = 0, 1, 2, and 






12 

is the photon propagator in the ^-gauge (the Coulomb gauge corresponding to ^ = 1 and the Feynman gauge 
corresponding to ^ = 0) with infrared cutoff at momenta of the order 2^ , h* < 0, then P{d^) and P^{dA) are 
the gaussian "measures" associated to the propagators S'o(k) and w^-'^ (p), respectively: 

^('^*)=]77 n n ^*k,.,p^*k,.,pcxp{--i^ Y. KAs^o-)v"K.} ^ (3-4) 

1 1 * 

pL'^'^dA)^-. n J] dRcA^_prfimA^pCxp{-— — ^ i^,p[z&^'''*(p)];Ji,,_p}, (3.5) 

where A/'ip, A/^ two normalization factors; in the second line, the product over p runs over the subset 'Dt j^ of 
Vp^L such that p > (here p > means that either po > 0, or po = and pi > 0, or po = Pi = ^^d p2 > 0) 
and p is in the support of x(|p|). Similarly, the * in the sum at exponent indicates that the summation runs 
over the momenta in the support of x(|p|). 

Remark. Note that, by the very definition of the cutoff function x (see the lines following Eq.(2.6)), the modes 
Ap.p associated to momenta close to ±(pj — Pp) and to their images over the dual lattice A* (i.e., to the points 
(0, ±(pp — p^) + TiiGi + ^26*2), with ni, ri2 G Z) are vanishing. In other words, the UV cutoff on the photon field is 
chosen so small that umklapp processes are suppressed. 

3. The interaction is 



X-\-Sj—Si,a,2 X,(T,1 



i=l,2,3 

Y J dx(ylo,x*+,,i*x,.,l + ^0,x+5i *+.,2*x,.,2) , (3.6) 



le 



/2 



where Jdx is a shorthand for /_«/2 ^^2:0 X^seA- "^^"^ interaction can be equivalently rewritten in momentum 
space, see Eq.(B.l). 

4. The external sources are 

s(*, A,^)^ j d^Y. E *?.xCj:x , (A, ^)= fd^Y (^x,.*x,. + K^Ka) , (3.7) 

a j = l,2,3 •' <T=ti 

where the sum over the index a runs over the choices a = K, CDW, AF, D, J, C, and Cf ^ are given by the same 
expressions Eqs.(2.15) after the replacement of the femionic and bosonic operators by the Grassmann and real 
fields ^ and A. The external fields A^o- a-re Grassmann (two-components) fields, while $",(. and J^^ are real 
fields, the first defined on the lattice, the second in the continuum with the same UV cutofi^ as the photon field 
(therefore, as discussed in the previous remark, the UV cutoff on the J field is chosen so small that the modes 
jfi,p with p sufficiently close to ±(pj^ — p^) and to their images over A* are vanishing). The source term can 
be equivalently rewritten in momentum space, see Eq.(B.l). 

The response functions introduced in Section 2 can be written as functional derivatives of the generating function, 

<(- - y) = A- A- .;i-^ .JLa ^'^'' ('J>, 0, 0)|,^„ (3.8) 



«,x" ],y 

that, remarkably, are independent of the choice of the gauge, ^ G [0, 1] A sketch of the proof of the functional integral 
representation of the observables of our theory, as well as a proof of the independence of Eq.(3.8) on the specific choice 
of the gauge, is given in Appendix A. Since the response functions are independent of ^, from now on we choose to 
evaluate them in the Feynman gauge, i.e., £, — 0. 

Similarly, the Schwinger functions in the ^-gauge at finite volume and finite temperature are defined as the /i* — ?■ — 00 
limit of 

'-'n,m;e,CT,p,;il^li • • • J^myii • • • ^ymj — ^,£1 ^>e„ p. j ^j x , n' ^ ' 

<^'^Xi,0-i,pi • • • t'/\x„,<T„,p„t'.^Pl,X„ + i • • • 'j'-^/J^.Xn + m A = ,y=0 
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Contrary to the response functions, the Schwinger functions (at least the way they are defined in Eq.(3.9)) are not 
gauge invariant and, therefore, they depend on the specific choice of the gauge. Nevertheless, we decide to compute 
them in the Feynman gauge, which is technically the simplest where to perform computations. Although the Schwinger 
functions in the Feynman gauge do not have an obvious Hamiltonian counterpart, we believe that they are a source 
of valuable information on the behavior of the system. In particular, as we will see, the 2-point functions and the 
vertex function, defined as 

[S2M (k)J ,,, = PL -p -^-— , ^§ 2;(^..) (P) = /3|5l I —7J-JJ ' 



mA^^P)],,-P'L-\SL\ ^ '-' , (3.10) 



will play a crucial role in the study of the flow of the effective couplings (and, therefore, of the response function 
themselves) . 

The independence of the gauge invariant observables on the specific choice of ^ is strictly related to the gauge 
invariance of the generating functional W^'^ with respect to U{1) gauge transformations. Namely, for all ^ G [0, 1] 
and h* < 0, we have: 

0=-?-W«^''*($,J + 5a,Ae*"") , (3.11) 

Oap a=0 

see Appendix A for a proof. By taking derivatives with respect to the external fields in Eq.(3.11), we can generates 
infinitely many identities between correlations, also known as Ward identities. In fact, Eq.(3.11) is equivalent to 

^^' dj ^Iis7\ ^ Vp,.ro(p) -^ + ^ ro(p)Ak^, , (3.12) 

with ro(p) = I „ -ipSi )• Taking, e.g., one derivative with respect to J or two derivatives with respect to A, we 
find: 

2 2 

EPA.4%^;k-)(P) = ' EPM>5|a;M(k,p) = e(ro(p)4:o*(k) - 5|f (k + p)ro(p)) , (3.13) 

(i=0 p=0 

which will be used below to deduce that the dressed photon mass is zero and the dressed electric charge is close to 
the bare one within O(e^). Note the crucial fact that Eqs.(3.11)-(3.13) are valid at finite volume, finite temperature 
and for any value of h* . 

4. RENORMALIZATION GROUP ANALYSIS 

In this Section we start the evaluation of the generating functional Eq.(3.1). In the following, for simplicity, we set 
the external fermionic and bosonic fields to 0, J = A = 0, and h* — — oo (dropping the h* index in the formulas). 
The effect of the external fields J and A has been discussed several times in the literature, see, e.g., [18, Appendix 
B] . The presence of a finite bosonic infrared cutoff on scale 2^ will be discussed at the beginning of the next section. 
Moreover, we set all the external fields $''"^ but the one coupled to the Kekule distortion to zero. The presence of 
the external field $(^) will be discussed in detail, for illustrative purposes. The effect of its addition is non trivial, 
particularly because of the new marginal terms it can generate. The effect of other external fields $("), a ^ K^ can 
be studied along the same lines and will be discussed in Section 7. 

As mentioned above, from now on we will work in the Feynman gauge, ^ = 0, in which case, the bosonic propagator 
simply reads 

M.)^.,:^.iMl^^l^j^y (4.1) 

The bosonic propagator w(p) is singular at p = 0, while the fermionic propagator is singular at the two Fermi points 
k = Pp. The first step of the RG analysis consists in rewriting both the fermionic and the bosonic propagators as 
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sums of two propagators, one supported close to the singularity (infrared propagator) and one in the complementary 
region (ultraviolet propagator), that is, 

S'o(k)=.9(^")(k)+.9(i)(k), tI;(p)=w;(^o)(p)+u-«(p). (4.2) 

where 

5(^")(k) = ^ x(|k - P^l)^o(k) =: Y: 5i^"'(k - P^) , w^^"\p) = X(IPI) / P^^ = ^ ■ (4.3) 

Note that in the first formula the support functions x(|k— pj|) and x(|k— p^|) have disjoint supports. Correspondingly, 
we rewrite the Gaussian measures as 

P(d*) = [ Yl P(d*£^"))]P(d*(i)), P^{dA) =P(dA(^"))P(dA(i)) , (4.4) 



where ^ir , 'i'^^\ ^f-"' and A^^^ have propagators given by glr , g'^^\ w'-"' and 'w'^^\ respectively. The fields 
^oT are called quasi-particle fields, and the index w = ± is called quasi-particle or valley index. 

Using this decomposition of the fermionic and bosonic fields, the generating functional >V($) = yv?'^°° ($, 0, 0) L „ 
at J = A = can be rewritten as 



„W(*) ^ / p(rf,5(<0))p(rf^(<0)) / p(-^^(l))p(-^^(l))gV(>I'<^°)+*<i),A(:^«)+A<i))+8(*(^")+*(i),A<^'')+A(i),*) ^ 



where the expression in the second line is obtained by an explicit integration of the UV degrees of freedom, which is 
very simple: in fact, the UV theory for the imaginary time variable in the presence of a fixed UV cutoff on the spatial 
variables is trivially convergent, see [15, Appendix C] or [48] for more technical details on this issue. The quantities Fq 
and S'*^-°''($) (normalized so that S'^-'^^O) — 0) are the contributions to the specific free energy and to the generating 
function of response functions, respectively, coming from the ultraviolet integration. V^''-* and B^^' (that are both 
normalized in such a way that V^°''(0, 0) = B^'^^{0, 0, <&) = 0) are the effective interaction and the effective source term, 
whose structure will be explicitly spelled out below. 

The integration of the IR effective theory is performed by using an iterative procedure, based on the following 
decomposition: 



5, 



,^°)(k') = ^£'')(k') , gW(k') = A(k')^o(k' + p^) , 



h<0 



hip) 



^I;(^°) (p) ^ E -^'^ (P) ' -^'^ (P) ^ ^ ' (4-^) 



/i<0 



where //i(p) = x(2~''|p|) — x(2"''^^|p|)- At each step we integrate the propagators g^'^^ and w'-'^\ h = 0, —1, —2, . . . , 
corresponding to degrees of freedom on momentum scale of order 2", 2~^, 2~^, . . ., the result of the integration defining 
the new effective interaction and source terms at scale h. At each step, we identify the marginal and relevant terms in 
these effective potentials and correspondingly define the effective coupling constants at scale h. Finally, after having 
inserted the quadratic fermionic relevant and marginal terms in the gaussian Grassmann integration (so defining a 
fiowing dressed fermionic propagator), we proceed to the next integration step. After the integration of the scales 
0, —1, . . . , —h + 1, we get (see below for an inductive proof) 

where P{d'i'ur ) and P{dA^-'^') have propagators 

^" ^^^ z4k')W(k'Mk'+p^) tko ; ' ^ ^p^ 2|p| ' ^^-^^ 

with X/i(p) = x(2~''|p|)and Z/j(k'), w^(k') the effective wave function renormalization and Fermi velocity at scale h, 
to be inductively defined below. Moreover, if Zh — Zh{0), the effective interaction and the effective source can be 
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written as sums over monomials of the fields (we recall that all the external fields <l''^"^ with a ^ K are set to zero, 
for notational simplicity): 



n,m>0: «— 1 i—1 



(4.9) 



n,m>0 

n+m>l 



271 Tn, jD 

n,?n>0: p>l i—1 i—1 i—1 

n+m>l 

where the integral sign is a shorthand for the sum over momenta and over the field labels, the underlined variables 
indicate a collection of variables (e.g., k = (k'j^, . . . , kj^)) and 5i^ enforces momentum conservation; note that, precisely 

because of momentum conservation, W'2„,m,p(k'j Pi q) explicitly depends on 2n + m + p — 1 variables rather than on 
2n + m +p (the "missing" momentum, which can be eliminated using the delta, can be chosen arbitrarily among the 

variables (k',p,q)). In Eq.(4.9) the kernels W^2„,m,p ^^^o depend on the choice of the field labels e,a,p,ui,ii, but we 
dropped these indices to avoid an overwhelming notation. 



A. Emergent relativistic theory 



In order to make the emergent relativistic structure of the theory apparent, it is convenient to rewrite the fermionic 
propagator in Eq.(4.8) as 



#'''(k') 



/Q r 2 \ k ,cj",(aJ k.' .ct.uj/ i-i^ 



1 



^h(k')ifco^Ji+^l;ft(k')fc'•^, 



■(l + i?,,,^(k')) 



where k' = (fco,fc'), |i?,j^^^(k')| < (const. )|k'| and 

r" = -1 , r^ 

with 



Cl 



1 

1 



(^2 = 



ICT2 



-i 
i 



iLoai 



CT3 = 



1 
-1 



(4.10) 



(4.11) 



(4.12) 



the standard Pauli matrices. We can introduce a "Dirac" 4-spinors, which makes the relation between the quasi- 
particle fields ^J and a theory of massless Dirac fermions more transparent: 



V^k'.iT V *k',o-,2,-' *k',cr,l,-' ^k',(T,l, + ' *k',cr, 



2, + / ' 



<f^' 






(4.13) 



The propagator of the tp, ip fields reads: 



^UlS<h)^(<h)^ _Xhiy) 



1 



where |i?;i(k')| < (const. )|k'| and 7^, /i = 0, 1,2, are euclidean gamma matrices: 



70 



1 
-1 



/ ia2 
^i^U 



72 



\iai 



ia\ 



(4.14) 



(4.15) 



satisfying the anticommutation relations: {7^,7,/} = —2(5^^. For what follows, it is also useful to define 73 — 

-io-3^ 
-10-3 



and the corresponding fifth gamma matrix: 



75 == 70717273 



1 
-1 



(4.16) 
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which anticommutes with all the other gamma matrices: {7^, 75} = 0, V/x = 0,1,2,3. 

Modulo the correction term R^, the propagator in Eq.(4.14) is the same as the one for euclidean massless Dirac 
fermions in 2 + 1 dimensions. The analysis is therefore very similar to the one performed in [18] for a system of 
interacting Dirac fermions coupled to with a massless gauge field. In particular, the scaling dimension of the kernels 
WQ„ m p of the effective potential is the same, see [18]: 

D = 3-2n-m-p, (4.17) 

where we use the convention that positive scaling dimensions correspond to relevant operators and viceversa. 

B. Lattice symmetries 

An important difference between the present case and the one studied in [18] is that here the propagator is not 
exactly equal to the Dirac one: on the contrary, it differs from it by the correction term proportional to i?;i(k'). Of 
course, this correction term is dimensionally negligible: therefore, it does not change the power counting. However, it 
violates some continuous relativistic symmetries used in [18] to exclude the presence of several relevant and marginal 
terms in the RG flow. One may fear that the lack of such symmetries might be responsible for the generation of new 
marginal or relevant terms, which are absent in the relativistic Dirac model. These potentially dangerous terms can 
be controlled through a careful analysis of the honeycomb lattice symmetries. In particular, it is proved in Appendix 
B that the effective interaction, the effective source and the gaussian integrations at all scales h < are separately 
invariant under the following symmetry transformations. Again, we spell out the symmetries only in the presence 
of the external field $^, the effect of the other external fields being discussed in Appendix B. In the following 
formulas, we drop the scale label h for notational simplicity; moreover, we think of ^j^, ^^ (^k' a u) ^^ being the the 

column (row) vector of components ^j^^, ^ „ ^ (^k' o- p w)' P ~ 1' 2' ^"^^ ^^ think of Ap as being the column vector of 
components A^^p, ji — 0,1, 2. 
(1) Spin flip : ^^'.^^c^ ^ ^I'-a^u,^ and A^^p, ^fp are left invariant; 

s-R" 



(2) Global t/(l): *k/ ^^^^ ^ e"'"''*^, ^^, with Q!,^ e R independent of k', and A^,p, $j^p are left invariant; 



(3) Spin SO (2): ( ^ ^'^f^''^" -^ e'^''^ ~ ^'■f''''" ) , with 9 independent of k', and A^^p,'l>f^ are left invariant; 



(4) Discrete spatial rotations: if Tk = (fco, e * 3 °'^k) and n_ = (1 — o-^j^ 

i p«(p'F+k')(^3-*i)n-Vl/^ \T/+ ^ \: 

^ ^ ^Tk',o-,cj ' ^k',a-,w ^ - 

-1 ^_ l^K . ^K 






Ap^T-'Arp, $-p ^ ^jVi^^p 

(5) Complex conjugation: if c is a generic constant appearing in the effective potentials or in the gaussian integra- 



tions: 



^ ' '- -I ^k',(T,Lj ' ^— k',CT, — Lj -I ^^P ' ^^— P ' 3,P 3, — V 



(6. a) Horizontal reflections : if i?ftk = (fco, — fci, fc2) and rul = 1, 7*^2 = 3, rui = 2, 

*k',.,.^'^i*H,k',p,., *k',.,. ^ *tk',.,c.'^i > Ap^R^A^^pC^^'^ $fp ^ df^^. ^^pe-P(^^-^^) ; 
(6.b) Vertical reflections : if i?t,k — (fco, fci, — fc2) and r^l = 1, r„2 = 3, r„3 = 2, 

^l,a,u -> *kk',.,-<. , Ap -^ R^Ar^p , ^f^p -^ $f„j,fl„p ; 

(7) Particle-hole : if Pk = (fco, -k), 

*k',.,c. ^ **Pk',a,-c. ' ^P -^ P^-pp ' ^L -^ ^f.-pp ; 

(8) Time-reversal : if /k = (—fco, fc). 

In the following subsection, the implications of these symmetries on the structure of the marginal and relevant terms 
are discussed. 
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C. Localization and the symmetry properties of the local terms 

In order to inductively prove Eq.(4.7), we write Vf*' == rvC*' + 7^V('') and B^''^ ^ CB^''^ + TZB'-''\ where the £ 
operator isolates the local terms, while TZ isolates the irrelevant terms: according to Eq.(4.17), we define 



W^u,m,p:cu,i^j (0,0,0), if 2n + m + p = 3 , 

^WZ'm,p;u,,^Ak',P, q) - <! [1 + (k', w) "%.,^) + p • ap + q • aq)W^i^!„,p.^,^,, (O, O, O) , if 2n + m + p = 2 , 



, otherwise 



(4.18) 
In the second line, (k ,uj_) ■ 5(k' (^) = '^iZiO^ii'^i) ' ^^(k'.wi), where (k',CLi) • 9(k',cj) is a shorthand for 

(k', u;) ■ a(k',c.) = kodk„ + nijTp + k') dp^^ + r!*(p^ + fc') dt^^ , (4.19) 

with dr, = -(—idk' +ujdk') , dp, = -(idi.' +ujdk') ■ 

In Appendix C it is proved that, thanks to the symmetry properties (l)-(8) listed in the previous subsection, the 
only non- vanishing local terms with 2n + m + p = 3 are either those with (2n, m,p) = (0, 2, 1) (i.e., terms of the form 
^^ AA) or the vertices, with 2n = 2 and m + p = 1 (i.e., terms of the form A^+^~ or ^^'i^'i^). The latter have 
the following explicit structure: 

'C<!o;(...).,(k',P) - *A,,,r^ , /:<o!'4,±.),,(k',P) == %r±^- , (4.20) 






where F, ■ := ' 



and the apex K added to the kernels with p 7^ is meant to remind the reader that the external field $ is of type K. 
The constants A^./i, Zj^y^, Zh are real and Ai^/j = \2,h- Note that the kernel with 2n = 2 and m = 1 with different 
omegas is zero simply because the photon field A, as well as the external field J, has an UV cutoff that makes the 
modes corresponding to momenta p close to ±(pj; — p^) and to their images over A* vanishing, see Remark after 
item 2 in Section 3. 

Moreover, the only non- vanishing local terms with 2n + to + p = 2 are either those with (2ri, m,p) — (2, 0, 0) (i.e., 
terms of the form \I'+\I'~) or those with {2n,m,p) = (0,2,0) (i.e., terms of the form AA). They have the following 
explicit structure (see Appendix C for a proof): 

^^^.°.o=(--)^^)=U.^(p^ + fc') ■ -0." j' ^^o,..o;(...)(P)=-../^'5.^, (4-21) 

where z^^h and Vfj__h are real and i^ui = i>2.h- 

D. The single-scale RG step: the inductive integration procedure 

The splitting into local and irrelevant terms is used in the inductive integration of the generating functional in the 
following way: we rewrite the integral in the r.h.s. of Eq.(4.7) as 

/"p(d*(-''^)P(dA(-''he^"''(^^**"'''"^*"''')+^"''(^^**"'''"^*"''''*) = 
^g/3i"th /"p(^^(</i))p(^^(<M)e^"''(v%:*<2''',A<S''))+e('')(v%:*<^'^',A<£''),<E.) ^ (-422) 



where V'^''' ~ V*-''-' — L^V^'^\ with L^V^'^^ the contribution to £V^'*'' that is quadratic in the fermionic fields (i.e., 
the one corresponding to the first term in Eq.(4.21)) and th is a normalization constant. Moreover, P{d'^^-'^^) has a 
propagator given by the same expression as Eq.(4.8) but for the fact that Zh(}^') and f/i(k') are replaced by Z/i_i(k') 
and w/i_i(k'), respectively, where 

Zh-i{\^) = Zft(k') + Z,,zo,ftX/.(k') , Zh-i{^)vh-i{\^) = Z;,(k')vft(k') + Z,,zi,^x/.(k') , (4.23) 
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and Zh = Zh{0), Vh — Vh{0). After this, defining Z^-i = Zh-i{0), we rescale the fermionic field by setting 

B'■''H^J%¥^^''\A'^^^''\'^) =: B''"'Hy^Z^¥-^''\ A'--^''\<P) . (4.24) 

By using Eqs.(4.20)-(4.21), the local part of the rescaled effective potential can be written as 

CV^'\^/z;r.^^^''\A(^'^^) = ^^Y.[Z'^-^^^^^J^^^^^^ (4-25) 

where e^,^ := -^^X^i^h and, if w;i_i := Vh-i{Q), 

'P ' Rl'^ ^ ^k'+p,<T,w^ w*k',CT,w ' Jp • ot2 Z_^ ^k'+p,<T,w^ "^k',cr,w • V^-'^^J 



W,(T 



k' w,(T.k' 



Note that, by using the notation defined in Eqs.(4.13),(4.15), the density and the current in Eq.(4.26) can be rewritten 
in the familiar relativistic form as 

°'P ■" /3i2 Z^ ^k'+p,o-70?'k',cr ' Jp •" 0-^2 / . ^k'+p,<T7yk'.g- ' l^-^'J 

(T,k' cr.k' 

Finally, by using Eq.(4.20) and the properties stated right before this equation, we find that the local part of the 
effective source term is given by 



k',p 

Note that, in contrast to what happens in a relativistic QFT, the effective source term Eq.(4.28) contains marginal 
terms that were not present in the original functional integral, i.e., the terms ^^ AA in the first line of Eq.(4.28). 
These potentially dangerous terms can be shown to be harmless by using the lattice symmetries (1) - (8), see Section 
5 B for a discussion of this point (nevertheless, let us anticipate that the reason why these terms do not create troubles 
is that they are "almost zero" , precisely because they vanish in the relativistic approximation; therefore, their naive 
dimensional bound can be improved and they can be shown to be effectively irrelevant). 
After the rescaling Eq.(4.24), we rewrite the r.h.s. of Eq. (4.22) as 



^-PLHn / P{^d^i<h-i))P[dA'^<'^-^^) / p(rfvl/('0)p(dA('^))eV"''(V^I^*'-'".^<-'")+^<'''(V^^*<-''^^'-'">*) , (4.29) 



where P((i*(^''"i)), P{dA'-^'''-^^) have propagators given by (4.8) with h replaced by /i- 1, while P{dA^^^), P((i*('*)) 
have propagators 

2|P| A-i A-i ifcoEO +zi;,j_i(k')fc' -r^^ 

where A(k') := Z;,_iA(k')/Z,.-i(k') and |i?^^^(k')| < (const.)|k'|. 

Remark. The single scale propagator can be decomposed as a sum of a Dirac-like propagator 5]j^^(k'), which is 
the propagator obtained by setting B!y^^ = in the second definition in Eq.(4.30), plus a rest, which has a better 

infrared behavior. We shall correspondingly write glj (k') = gjy^{kl) + ri, (k'). This decomposition will be useful 
in the following, as already anticipated by the comment after Eq.(4.28). 

At this point, we can finally integrate the fields on scale h and, defining 

= g-/3L^(i=^;.+tO + S<^'''(*) /"p(^v[r(''))P(dA(''))e*"''(^^^**"''''^'"'''^+''"''(v^^^**"'''''^*"''''*) , 
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our inductive assumption Eq.(4.7) is reproduced at scale h — 1. Note that Eq.(4.31) can be thought as a recursive 
definition for the effective potential. The integration in Eq.(4.31) is performed by expanding in series the exponential in 
the r.h.s. and by integrating term by term with respect to the gaussian integration P{dil)'^^'>)P{dA'^'^^) . This procedure 
gives rise to an expansion for the effective interaction and source terms in terms of the renormalized parameters 
{^tj.,k,Vtj.,k,Zk-iTVk-i,Zj^^,\Jj^}h<k<o, which can be conveniently represented as a sum over Gallavotti-Nicolo 
(GN) trees [14]; the value of each GN tree can be thought of as a sum over connected labelled Feynman diagrams, 
i.e., every GN tree represents a set of Feynman diagrams characterized by the same hierarchical structure of the scale 
labels associated to the propagators, see [18, Section 2.2] for a thorough discussion of this expansion. 

We will call {e^.fe, i'^ji,^k\k<o the effective couplings or running coupling constants. The constants e^^/i play the role of 
effective charges, while Vf^^h play the role of effective photon masses. It will be shown below that the effective charges 
stay constant under the RG flow (more precisely, e^^^h sue essentially independent of fj, and h) and that the effective 
photon mass is vanishing (more precisely, I'p.h is small uniformly in h, on the "right scale"). 

Remark. Note the unusual dependence of the effective charges on the index /x: the global symmetries (l)-(8) 
discussed above ensure that eijj — e2.h but they do not a priori guarantee that egjj = eijj. This situation is in 
striking contrast with what happens in QFT, where Lorentz invariance guarantees such a property to be valid at all 
scales. However, in the next section we will show that, thanks to lattice Wis, eoji and ei jj, even if not exactly equal 
to each other at all scales, admit the same limit as /i — > — oo; namely, eo,-cx) = ei,_co = e + O(e^). 

The expansion in GN trees and labelled Feynman diagrams allows us to obtain the following inductive estimate 
on the kernels W2n,m,p in Eq.(4.9). Let Sh — max/i<fc<o{|e^jt|, |f^jt|} be small enough. If Zk/Zk-i < e'-^^'' and 
C~^ < Vk-i < 1, for all h < k < and a suitable constant C > 0, then the N-th order contribution to W2n,m,p in the 
effective couplings (to be denoted by W^2n'.m'p) admits the following bound (the "iV! bound"): 

ll<:l^!pll<(const.)^£t(Y)!, (4.32) 

where ||VF2n'mpll •— 2^'''^^^^"^™^^) sup |W2Ti'm p(k'' PiQ)! ^nd the sup is performed with respect to the momenta 
and the field labels. The proof of Eq.(4.32) can be found in [18, Section 2.4]. The basic ingredient in the proof is a 
dimensional estimate of the kernels, which follows from the bounds: 

\w^''Hp)\ < (const.)2-'^ , ^ Y. I*^'^(P)I ^ (const.)22'' , 



p 



ll£'^^ 



(k')ll < (const.)2-'' , -L J2 \\eHk')\\ < (const.)22'^ , (4.33) 



that is, every propagator on scale h is associated to a factor 2^'' and every loop integral on scale h is associated 
to a factor 2^^^. The estimate Eq.(4.32) follows from: (i) counting the number of propagators and loop integrals on 
each scale h for every given labelled Feynman diagram; (ii) realizing that the corresponding dimensional estimate is 
uniform in the Feynman diagram, for all the diagrams associated to the same GN tree; (iii) performing the sum over 
scale labels for every fixed GN tree; (iv) counting the number of Feynman diagrams associated to each GN tree and 
the total number of GN trees contributing to order N in renormalized perturbation theory. See [18, Proof of Theorem 
2.1] for more details. 

The bound Eq.(4.32) tells us that the A^-th order contribution to the effective potential is finite in norm, uniformly 
in h. If the effective couplings remain small in the infrared, informations obtained from our renormalized expansion by 
lowest order truncations are reliable at weak coupling. The importance of having an expansion with finite coefficients 
should not be underestimated. The naive perturbative expansion in e the fine structure constant is plagued by 
logarithmic infrared divergences and higher orders are more and more divergent. More precisely, one can find classes 
of diagrams of order A^ in e contributing to the effective potential on scale h whose size grows like 0(|/i|"). Therefore, 
finite order truncations of the naive perturbation theory do not give a priori any reliable information on the IR 
behavior of the theory, not even at weak coupling. 

Regarding the combinatorial factor in the r.h.s. of Eq.(4.32), we note that the {N/2)\ dependence is compatible 
with Borel summability of the theory. However, summability does not follow from our bounds. Constructive estimates 
on large fields for the bosonic sector (in the spirit of, e.g., [6]) combined with determinant estimates for the fermionic 
sector (in the spirit of, e.g., [15]) may allow a full non-perturbative construction of the theory. However, this goes 
beyond the scope of this paper. 

Let us conclude this section by adding a comment, which will be useful for the study of the flow of the effective 
parameters discussed in the next sections. Consider the splitting guj (k') — gj-,|^(k') + ri, (k') mentioned in the 
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remark after Eq.(4.30). As mentioned there, the rest rij ' is better behaved in the infrared than the relativistic 
propagator gjj ^ . More precisely, 

||r('')(k')|| < const. , ^^||rW(k')|| < (const.)23'^ , (4.34) 

P^ k' 

which should be compared with the second line of Eq.(4.33). It is apparent that ri, is associated to a dimensional 
gain 2'* as compared to the leading term gjj^ik.'). This implies that if we decompose VFg^'^ „ as 

where W^^\^' is obtained from W^^mp ^y replacing all the propagators gij '(k') by (f^'^i^) (and by neglecting 
the contributions coming from the UV propagators on scale ft. = 1), then W^^mn i*^ dimensionally negligible in the 
IR as compared to the "relativistic" contribution W^^^p ! i-^-' ^2nmp admits a bound similar to Eq.(4.32), with 
an extra factor (dimensional gain) proportional to 2 , with < ^ < 1. This follows from the improved dimensional 
estimate on the propagator ri,' , Eq.(4.34), and from the fact that "long GN trees are exponentially depressed", i.e., 
the property referred to as "short memory property", see [18, Section 2.4]. 



5. WARD IDENTITIES AND THE FLOW OF THE RENORMALIZED PARAMETERS 

We have seen that the effective potentials (and, similarly, the correlation functions) can be written as series in the 
effective charges e^.^ and the effective masses i>fj,^h with bounded coefficients at all orders, uniformly in the infrared 
cut-off, provided that the ratios Zh/Zh-i remain close to 1 and that effective Fermi velocity remains bounded away 
from zero along the RG flow. Of course, such expansions are useful only if the running coupling constants remain 
small for all values of /i, a fact that we are going to prove to be true, thanks to exact lattice Wis. In this section 
we flrst study the flow of the running coupling constants {e^^/j, v^^hjhKa and next the one of the other renormalized 
parameters. 



A. The flow of the electric charge and of the photon mass 

The key idea is to get informations on the running coupling constants {e^^/i, i^^^ft};i<o by using the Wis Eq.(3.13) 
for the sequence of reference models W"''' ($, J, A), where the scale of the bosonic IR cutoff h* is thought of as a 
parameter. For each choice of the IR cutoff /i*, the generating functional W^'^ ($, J, A) is computed by a multiscale 
integration procedure similar to the one described in the previous section, with the important difference that after the 
integration of the scale h* we are left with a purely fermionic theory, which is super-renormalizable: in fact, setting 
TO = in the formula Eq.(4.17) for the scaling dimension of the kernels of the effective potentials, we see that the 
scaling dimension of the reference model below the cutoff h* is D ^ 3 — 2n—p. In particular, the kernels of the effective 
interaction (i.e., those with p = 0) are always negative once the two-legged subdiagrams have been renormalized; as 
shown in [15, 20, 21], the effect of the integration of the scales < h* is just to renormalize by a small flnite amount 

the effective parameters Z/j. , w/j. , Z^ ^, , X^ ,j. . Let us denote by {e' f^ , i' i^}h<o the running coupling constants of 
the reference model with infrared cut-off on scale h* . Of course, if ft > ft*, 

i^fiJi ' '^fi]h }h*<h<Q = {e^,h, l^f^,h}h'-<h<Q , (5.1) 

where the constants in the r.h.s. are those of the model without IR cutoff, i.e., e^ji :— e ^ and ly^ji :~ v j^ . On 
the other hand, as proved in [18, Appendix B], the two- and three-points correlation functions of the reference model 
are proportional to the inverse wave function renormalization and to the effective charges, i.e., if k' and p are such 
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that |k'| = 2'^ and |p| < 2'' , then 



4°:o'(Pf + k') = ^|M(i + B^,,.(k')) , (5.2) 



Zh'- 



1 



4°f;o(PF + k', P) ^ ^ ^"' ^^''^ + P^ (eo,,.r° + ei?:^,,, (k', p))gi''^)(k') , (5.3) 

4°f.;(p^ +k',p) = ^vy,,^^^-^^L_^±^Un, rUeSL..*(k',p))e*Hk') , ? e {1,2} , (5.4) 

where the correction terms i?tj^^.(k') and i3^^.(k',p) are of order e\, (recall that Sh = niax{|e^^fc|, |j^^,fc|};i<fc<o), 
uniformly in /i*, for all |k'| ~ 2^ and |p| < 2'* . For later use, let us note that B(^^/j*(k') is difFerentiable in k' and 
its derivatives computed at |k'| — 2^ are dimensionally bounded by (const. )2^'' e^.. 

Thanks to Eqs.(5.2)-(5.4), we see that informations on the mutual relations between the correlation functions of 
the reference model with cutoff h* (which are provided by the Wis) imply relations between the effective charges on 
scale h. Regarding the photon mass, an equation similar in spirit to Eqs.(5.2)-(5.4) is valid, namely, if |p| < 2 , 



4,2.(,,.)(P) - 2'^ ['^,M'S,. + <nP)j , (5.5) 

where the correction term i3^r(p) is of order e\,, uniformly in h* , for all |p| < 2^ . The proof of Eq.(5.5) can be 
worked out along the same lines of [18, Appendix B] and is based on the following remarks: for all scales h > h* , hy 
construction, the external field J appears in the effective potential in the combination A+J, see Eq.(3.1). As discussed 

in the previous section, the corresponding kernel, M^g 2 of*' t')^^-' ^^ equal to Wq 2 o(u iy)(P) ~ 2'* ( t'p ji(5p^ + i?^''(p) j, 
with i?|^''(p) a bounded correction, of second or higher order in the effective coupling constants. Therefore, after the 
integration of all the scales h > h* , we are left with a kernel J J equal to the r.h.s. of Eq.(5.5), with a slightly different 
correction term B^, (p) replacing i?,'^^ (p) . From that scale on, we are left with the super- renormalizable theory studied 
in [15, 21], in which the only marginal interactions are the J^+^^ terms, whose coefficient is renormalized by a finite 

amount under the RG flow from h* to — cx), see [21]. Therefore, the dominant correction terms to S^.'^./ ^\ip) coming 
from the integration of the scales below h* are obtained by contracting two effective vertices J^+^~ on scale h < h* , 
and then summing over h. The resulting contribution is dimensionally bounded as X]/i<ft* 0{2'^ef^,) = 0{2^ £^*)- 
Higher order corrections are bounded in a similar way, using the hierarchical structure of the GN trees, see also [21]. 
Now, combining Eq.(5.5) with the first WI in Eq.(3.13) computed at p = pn^ (with u^ the unit vector in direction 
/i and p = jpj < 2^ ), we find: 

v^^h' = - limB^^^ipu^) - 0{el) , (5.6) 

which means that there is no spontaneous generation of the photon mass (i.e., the photon field remains unscreened). 
If read in naive (non-renormalized) perturbation theory, the above identity is equivalent to an infinite sequence of 
cancellations taking place at all orders among the graphs contributing to the photon mass. At lowest order, the 
cancellation takes place between the two graphs in Fig. 3, as discussed in Appendix E 1. 

FIG. 3. The lowest order contributions to the "mass" of the field A^ with p = 1, 2, which cancel out exactly when computed 
at transferred momentum p = 0, as proved in Appendix E 1. 

A similar argument can be applied to control the fiow of the effective charge. In fact, note that by computing the 
second WI in Eq.(3.13) at (k', p) = {2^ Up,pUp) and by taking the limit p ^ 0, we find: 

4£(PF + 2'''u^,0) = -ea^4"f (p^ + 2''*u^). (5.7) 
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By plugging Eqs.(5.2)-(5.4) into Eq.(5.7), we get 

eo,h' - efl + 5° ..(Pf + 2'^*Uo, 0) - 2'**aoB^,^.(p^ + 2''*Uo)) =: e(l + Aom-) = e{l + 0{el)) , (5.8) 



eiM' - e(l + Si_;,.(p'^ + 2'**ui, 0)ri - 2""' diB^,h'{PF + ^'"'^i)) -: e(l + ^i,;,.) - e(l + 0(e^)) , (5.9) 

which tell us that the effective charges e^^h* remain close to the unperturbed value e^^ = e at all scales h* < and 
at all orders in renormalized perturbation theory. If read in naive perturbation theory, Eqs.(5.8)-(5.9) are equivalent 
to infinitely many cancellations taking place at all orders among the logarithmically divergent graphs contributing 
to the dressing of the electric charge. At lowest order, the cancellation takes place between the graphs in Fig. 4, as 
discussed in Appendix E2. 





[.S^.^] 



= 



FIG. 4. The lowest order contributions to the dressing of the electric charge, which cancel out exactly when computed at the 
Fermi points k = p'^ and at transferred momentum p = 0, as proved in Appendix E2. 

The corrections A^_h-' in Eqs.(5.8)-(5.9) are given by sums over GN trees of order two or higher in the effective 
couplings; their iV-th order is bounded by (const. )^(e/i.)^(A^/2)! Moreover, \A^^h* ~ ^/^,-co| = O{e'^^{v^oo ~ Vh)) + 
0(e?.o^2''/^), see [18, Section 4 and Appendix B] for a proof of these facts. Note that a priori Aq.-od t^ ^i,-oo; 
however, the approximate Lorcntz invariance of the theory combined with the flow equation for the Fermi velocity 
will allow us to show that j4o.-oo ~ j4i._oo, see the next subsection for a discussion of this point. 

Remark. The discussion in this section is very similar to the corresponding one for a model of Dirac fermions 
with electromagnetic interactions in the continuum in the presence of a fixed rotationally invariant UV cutoff, see 
[18, Section 3]. The key difference is that the analogues of Eqs.(5.8)-(5.9) in [18, Section 3] have an extra term in 
the l.h.s., i.e., they read e^,/j.(l - a^) = e(l + A^,^.), with \A^^h' - -4^,-oo| = 0{e'i^{v-^ - Vh)) + 0{e^_^2''/^) 
and |Ao,_oo — ^i,-oo| < (const. )e^; see [18, Eqs.(4.13)-(4.14)]. The extra constants a^ in the l.h.s. are due to the 
corrections to the Wis produced by the momentum cutoff used in [18], which explicitly breaks gauge invariance; 
as shown in [18], ao 7^ Q^i, the difference being of second order in the electric charge, which implies that Lorentz 
invariance is not recovered asymptotically in the IR. This has to be contrasted with the case studied here, where exact 
lattice gauge invariance is preserved by the RG flow, so that no correction terms appear in the l.h.s. of Eqs.(5.8)-(5.9) 
and Lorentz invariance spontaneously emerge in the deep IR, as proved in the next subsection. 

B. The flow of the effective parameters 

We are now left with studying the evolution of the effective parameters Vh, Zh^ Z^ ^, \^ ^ under the RG flow. If 
we use the analogue of the decomposition Eq.(4.35), we can write the flow equation for these parameters as: 

"^h^^l + Pl + rl, ^ = l^Pl + rl, (5.10) 

Vh Zh 



^K.h 



Ph^'''h ' ^j,p,,h-l ^ ^j,fi,h + ''^j,fi,h ^ (5-11) 



where /3^ are the relativistic part of the beta function that, by deflnition, are obtained by replacing all the propagators 
gli (k') contributing to the r.h.s. of the flow equation by their relativistic part gjj^(k'), while r* are the rests, which 
are smaller by a factor 2^'', 9 € (0, 1), as compared to the corresponding dominant terms (the reason is the same as the 
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.# 



one sketched after Eq.(4.35)). Note that both /3* and r* are functions of the whole sequence of effective parameters 
{ef,,k,'^ii,k,Zk-i,Vk-i,Z^i^,Xf^^j^}h<k<o that are bounded respectively by 0(£^) and 0(2'"'e^), provided that the 

ratios Zk/Zk-i are close to 1 and that Vk are bounded away from zero, for all h < k < 0. 

In the second equation in Eq.(5.11), we used the fact that the dominant part /3^ ^^ is exactly zero for all choices of 

j, ^, h, by a parity argument (inspection of perturbation theory immediately shows that all the contributions to (3; 

are given by integrals of odd functions of k' over a domain that is invariant under k' — > — k'). The bound on 

immediately shows that Xf ^ is uniformly bounded by 0{e1) in the IR. 






Regarding the flow of Vh, Zh, we note that modulo the correction terms r^^ and rj 



. f^ , they are the same as those for 
Dirac fermions in the continuum, derived and written down in [18]. In particular, using [18, Eq.(3.14)], we can write: 



Vh-l 



= 1 



log 2 

47r2 



:e\l - Vh){l + ^'J + -e(l + B'^){e^,h - e^^h) 



+ rl 



(5.12) 



where AJ^ is a sum of contributions that are finite at all orders in the effective couplings, which are either of order 
two or more in the effective charges, or vanishing at w^ = 1; similarly, B'^ is a sum of contributions that are finite at 
all orders in the effective couplings, which are of order two or more in the effective charges. The proof of Eq.(5.12) is 
based on the remark that fih_y would be vanishing if Vh '■= 1 and eo,/i '■— ei^h, by Lorentz invariance (see [18, Section 
3.3] for more details); the numerical coefficients follow from an explicit computation, see [18, Appendix C]. From 
(5.12) it is apparent that v^ tends as ft, — > — oo to a limit value 



w-co = 1 + 7^(eo,-oo - ei._oo)(l + C"_ao) 
be 



(5.13) 



with Ci^oQ a sum of contributions that are finite at all orders in the effective couplings, which are of order two or more 
in the effective charges. The fixed point (5.13) is found simply by requiring that in the limit h — > — oo the argument 
of the square brackets in (5.12) vanishes. 

Consider now the identities Eqs.(5.8)-(5.9). In a Lorentz invariant theory (i.e., in a theory where all the propagators 
guj (k') are replaced by their Dirac approximations g\j i;(k'), the Fermi velocity is equal to the speed of light, Vh '■— 1, 
and the charges e^^h are /x-independent, cg.h — ei.^), we would get that Aqji = Aiii. Therefore, using the GN trees 
representation for Ao,h, Aij^, their approximate Lorentz symmetry and the short memory property (see [18, Section 
2.4]), we find: 

eo,/i_i - ei,h-i = Ei,h{eo,h - ei,h) + £^2,h(l - Vh) + E^.h , (5-14) 

where Ei^h and E2,h are 0(£^) and E3,h = 0(6^2'"'), for 9 e (0, 1). If we combine Eq.(5.14) with Eq.(5.12) we get 

eo,-oo — ei,-oo = -E'-co(eo,-oo — ei,-oo) , (5.15) 

with E_oa = 0{e'^^), which implies the spontaneous emergence of Lorentz invariance in the deep IR, i.e., 

eo,-oo = ei,_oo and V-oo = 1 • (5.16) 

Using these informations into Eq.(5.12), we see that the approach of Vh to the speed of light is anomalous, i.e., 



2e' 



1-VhO^ A{v) 2'^\ fj^ +0{ 



(5.17) 



where A{v) is a function of e and v that vanishes linearly at u = 1, i.e., for v close to 1, A{v) = (1 — u)(l + 0(1 — 
v) + 0{e^)). The "~" in the first equation means that the ratio of the two sides tends to 1 as ft ^ — oo. 

Regarding the flow of the other renormalization parameters, a second order computation of the dominant contribu- 
tion to the beta function (see [18, Appendix C] and Appendix D 1), the fact that e^^h are close to e_co = e(l + 0(6^)), 
asymptotically as ft — >■ — oo, together with Eq.(5.17), shows that 



/?!- 
Pt- 



127r2 
127r2 



log2 + Oie^il - v)2''"'') + 0{e^) , 
■ log 2 + 0(e2(l - v)2^''''') + 0{e^) 



Ph 



Ze^ 

47r2 



log2 + 0{e\l - w)2^'= '^) + 0(e*) 



(5.18) 
(5.19) 
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for some c > 0. Eqs.(5.18)-(5.19) imply that 

Zn ^ B" {v) 2-''" , Z+ ,^ ~ B+ (v) 2-"^'^ , Z^,^ ~ B' (v) 2'^^'^ , (5.20) 

with 

and B'^{v) = 1 + 0(1 — f ) + O(e^), which concludes the study of the flow of the renormalized parameters. 

6. THE KEKULE RESPONSE FUNCTION 

In this section we explicitly compute the Kekulc response function (the other responses can be obtained in a similar 
way and will be discussed in the next section), 



o2 

ii-lf' (x) == lim lim lim __ -_>V"''^' ($, 0, 0) . (6.1) 



$=0 



The iterative construction of the generating functional described in Section 4 induces an explicit representation for the 
Kekule response function in terms of GN trees, completely analogous to the one described in [21, Proof of Proposition 
1]. In particular, using the analogue of the decomposition in Eq.(4.35), we can rewrite 

i?!f(x)=<)'"(x)+i?(f)(x), (6.2) 



where R\^ (x) is obtained by replacing the lattice by the Dirac propagators in the expansion for R\: and R\, (x) 
is the rest. Using the same strategy leading to the bounds Eq.(4.32) and, more specifically, to [21, Eqs.(2.81)-(2.84)], 
we find that the A'^-th order contribution in renormalized perturbation theory to Rl, (x) is bounded, for all M > 0, 

by 








2)' ,^ ^^V Z„ y " ^ l + (2''|x|)^ 

for a suitable constant Cm- The factors 2'*2'^'*[l + (2''|x|)*^] represent the dimensional bound on all the labelled 
Fcynman diagrams corresponding to the same GN tree, where: (i) h is the lowest among the scales of the propagators 
in the diagram; (ii) h is the lowest among the scales of the propagators in a path connecting the two special vertices of 
type <I>-^v]/+v]/-; (iii) 2^* = 2''('^~2n-m-p) | jg ^j^g scaling dimension of the graph; (iv) 2^'^ is the dimensional gain 

coming from the fact that the locations x and of the external fields (^^ are fixed rather than integrated over the 
whole space-time domain (as it is the case for the Feynman diagrams contributing to the thermodynamic functions, 

where all the space-time labels of the vertices are integrated over the whole space); (v) [l -I- (2''|x|)*^] is the decay 
factor coming from the propagators on a path connecting the two special vertices of type $^4*+^^. Now, picking 
M — b, exchanging the order of summations over h and h, and summing over h gives (recall the asymptotic relations 
Eq.(5.20)) 



l-Ry'' '■ (x)| < (const.) |e| (^yj! ^ — 



- ■ (2''|x|)E 



|i?;^^(^)'^(x)| < (const.)^|e|~f^)! , (6.4) 



where we used the fact that rj^ > rj^. The correction term R^- (x) admits a similar bound, with a dimensional 
gain factor that implies a faster decay in real space. Using the symmetries in Appendix B together with relativistic 
invariance (see symmetry [18, (7)]) and proceeding as in Appendix C, we find that the symmetry structure of i?j,- (x) 

is: 

Rij (x) = a]^ 'e -^ — ■ — .._c(K) + faster decaying terms, (6.5) 
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where the faster decaying correction terms come from: (i) the irrelevant terms that, scale by scale, produce corrections 
smaller by a factor 0(2'*) as compared to the dominant terms; (ii) the Lorentz-symmetry breaking terms, i.e., the 
terms proportional to 1 — ii/i that conic from a rewriting of the effective Fermi velocity Vh in the definition of gD,uj as 
u/j = 1 — (1 — f/j); these produce corrections that, scale by scale, are smaller by a factor 0{e^2^'^°'^^^-^'^ '') as compared 
to the dominant terms. In Eq.(6.5), S,'^-^^ = "2,1]]^ — 2?7 = |^ + 0{e^) and a)f is a suitable constant, bounded 
in absolute value by (const. )'^( ^ ) ■ The 0-th order constant, which gives the dominant contribution to the Kekulc 
response function, is given by the value of the graph in Fig. 5. 
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FIG. 5. Leading contribution to the Kekule response function; a sum over ft < is understood. 
More explicitly. 
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where hWk — max{/i, k}. Recalling that F^ 



> + 0(e |x| ^'^ ) + faster decaying terms , (6.6) 



e-'"^"^(Ti = e*P'^(^j-*i)cri and using the fact that Tr{FJ;;CTiF^„CTi} = 2 



-iu:Oj 



for all V e {0, 1, 2} and that Tr{F^(TiF'^^cri} = for /i ^ v, we can rewrite 

(ik I' dp 



MK),Di 



i?-'-'^(x) = 4 E E^ 



-ip';^(x-g,+g,) y^KMVk) 



h.k<0 w 



ZhZk J 2n\B\J 2n\B\ 



//.(k)/fe(p)e 



^(k— p)x 



k P 

|kP|p|2 ' 



(6.7) 



modulo corrections 0(e2|x| ^+'> ) or decaying faster than |x| ^+'> at infinity. Using Eq.(5.20), we can replace 
^^z'uzf ^y (§Tr^)^2-2';K(''Vfc)2»?/i+»?fc^ modulo faster decaying corrections. Finally, rewriting, e.g., 2"^^ = |x|-''[l + 
((2''|x|)'' — l)], it is easy to realize that the contributions to the response function coming from the terms like 
|x|-''[(2'*|x|)'' - l] are 0(e2|x|-4+«"^'), see [48, Section 4.3.2]. Therefore, using the identity E/i<o A(k) = x(k), 



4'"'"w = 4(|o|^)%i'^'''^-''^E^-"^^^^-''^''^ 



2^|Sr^ ' |k|2 / 27:\B 



'^^ X(p)e-^P^^ , (6.8) 



modulo the aforementioned corrections. Finally, using the fact that / 2^m\X0^)G^^^j^ — ^^|^T^+ faster decaying 



terms, wc can rewrite: 



R. 



a;{K),D,, _ 



■ij 



w-^ 



27 / B- (w) \ 2 cos [pp (x - 6^ + S 



87r2VB0(w) 



.|4-f(») 



|kp ~ 877 |3c|^ 



+ 0(e |x| ^^ )+ faster decaying terms . (6.9) 



Using the fact that B*{v) = 1 + 0(1 - w) + 0{e^) we finally get Eq.(2.20). 

As we already remarked in the introduction, the effect of the interaction is that of enhancing the response function 



uiK), 



of the Kekule distortion, but still preserving its integrability, so that R\, (p) is finite. However, dpR\, (p) is singular 
in p = p"^. In fact, the Fourier transform of Eq.(2.20) gives: 

pj(^.-^.)|p_p+|i-e<'^' +e"'P^(^'"^^-)|p-p^|i"^"''] +more regular terms , (6.10) 

which implies that d^R\'^\p) diverges like |p — p"^]^^ at the Fermi points, as claimed in Section 2C. 



i?(f(p) = (const.) 
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7. OTHER RESPONSE FUNCTIONS 



An analysis analogous to the one for the Kckulc response function can be worked out for the other responses and is 
sketched below. Regarding the CDW and AF responses, the symmetries of the model imply that the relevant terms 
of the form $'^^^A or ^^^A are vanishing, while the marginal terms of the form ([)<^^^vlf+vl/- and ^-^^vp+vl/- have 
the following structure: 



LB^J^L.i^Z^Ti^^-^'lAC^''^: 



■'CDW^ 

' y 

k ,p 



$) = 



\7+ S,CDW^{<h) + 



■^ K.a.u. 



Zkh'^ 



CDW 



* 



(<h)H 



xe'-'^^-.wr- 



(7.1) 



£^i'^](v/z;:r7vi/(<'o,A(^M^co) 



1 



P^L^\Sl\ 



k',p 



■J k',(T,a; 



'K.h 



$ 



AF 



* 



(<h) + 



\-p,a.uj L 



-3]^!^'^'- 



k'.CT.- 



(7.2) 



Similar expressions are valid for the other external fields $". Note that besides the marginal terms spelled out in the 
previous equation, there may also be marginal terms of the form $Aj4, $$A or p^$y4; the fiow of the corresponding 
coupling constants can be easily controlled along the same lines followed to prove the boundedness of the flow of \^ ^ 
in Section 5 B. 

As seen in the previous section, the flow of the renormalization constants ■^c/jvi'h ^^^ ^afh control the large 
distance decay of the CDW and AF response functions. In Appendix D 2 it is shown that 



Z. 



CDWJi-l 

7+ 
^CDW,h 



l + £2l°g2 + 0(e4) + 



^CDW.h-l 



z. 



CDW,h 



Vl-K- 



■log2 + 0(e4) 



(7.3) 



which implies that 



'CDW.h 



B. 



± 

CDW 



{v)2- 



^VCDW 



(7.4) 



with 



3e2 



VCDW 
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0{e^ 



(7.5) 



and B^jjy^^{v) = 1 + 0(1 — v) + 0{e^). The exponents for the AF are exactly the same, simply because the kernels 
of the marginal source terms ^-^^^Jz+x]/" are the same as the corresponding terms of the form (f'^-^^vp+vl/-^ see 
Eqs.(7.1)-(7.2). Therefore, proceeding as in the previous section, we get, for a ~ CDW,AF, 



«!;'w^^'^^-'"'^^ 



8tt^\B°(v)J |x|4-«<"' 



1 



^(CDW) ^ ^(AF) ^ 2{r]^ 



0(e^|x 



4e' 



2|^|-4+C* 



faster decaying terms 



CDW -'l) = ^ + 0{t 



(7.6) 
(7.7) 



As for the Kekule response function, the anomalous decay in Eq.(7.6) implies the presence of a singularity in the first 
derivative of the corresponding Fourier transforms, in analogy with Eq.(6.10). 

It is now clear how to extend such analysis to all other responses. The large distance asymptotic behavior depends 
on the specific values of the critical exponents, which are determined at leading order by a second order computation 
along the lines of Appendix D: one first identifies the structure of the local part of the effective source terms, in analogy 
with Eqs.(4.28),(7.1),(7.2); for each such term two renormalization constants appear, corresponding to processes with 
transferred momentum equal to or to p"^ — p^"^; the flow of each of them is given by a diagram like the one in Fig. 6 
and can be computed as in Appendix D. The resulting values of the critical exponents ^(c'laniei) _ 2{richannei — rj) ^^ 
second order in the electric charge are summarized in Table I. 

The interaction removes the degeneracy of the critical exponents and, therefore, we can identify the excitations 
whose response functions decay slowest at inflnity, which correspond to the order parameters in putative strong 
coupling broken phases. In particular: 
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channel 


critical exponent 


"■t,abuj,<y + bi,aO.'^,T 





a,uj,ab-ui,cr + b^^^a-ui,a 


4eV(37r2) 


i0.t,abui,a ~ ib'^,aa.Lj,r, 





ia-t,ab-L^,cr - ifei_CT«-cJ,<T 





0-t,,a'^i^,ty + bt,^abu>,<T 





0'uj,fT^ — ^,(7 ~r *^uj,a*^ — uj,(7 





at,aai^,a - bt,abui,a 


4eV(37r2) 


^uj,a^~^,(7 t)^^ fjO — Lt;,(T 






channel 


critical exponent 


O'iat^aO'^.tT + &i,CT&t^,<T) 





o"(ai,<Ta-cj,<T + b^^^b-^^a) 





O'iat^aO'Ui.a - &i,CT&t^,<T) 


4eV(3^2) 


o"(ai,<Ta-tj,<T - &i,CT&-t^,<T) 





^(ai,<Tat^,<T - &i,CT&t^,<T) 


4eV(3^^) 


ia(a+ ,a+ „^ + 6+ ,?)+_,) + c.c. 


-eV(3^2) 


iauj{a+^„a+_^ + b+^^b+_^) + c.c. 


"cV(3^^) 


<<T&-^,<T - <^<^ + c.c. 


-cV(3-^) 


jcr(a+ ^6l^^__^ - bi,<Talt^,-a) + c.c. 


-cV(3-^) 



TABLE I. Lowest order contributions to the anomalous exponents. Summation over the valley and spin indices is understood. 

1. the dominant exponents correspond to the second and seventh channel in the left table and to the third and 
fifth channel in the right table. The first three are the K, CDW and AF local order parameters discussed 
in the introduction. The fifth channel in the right table was introduced in [27]. Microscopically, this order 
parameter may be understood as a specific pattern of circulating currents and its enhancement in the presence 
of interactions is in agreement with the expectation that a time-reversal broken fixed point should emerge in 
the strong coupling regime [30]. 

2. The last four channels in the right table are Majorana-typc masses that correspond to intcr-nodc and intra-node 
(uniform and non-uniform) Cooper pairings [37]. Their critical exponents are all negative, which means that 
superconducting order is unfavored at intermediate to strong coupling. There are other possible Cooper pairings 
besides those explicitly reported in the table and it turns out that all their exponents are either equal to or 
smaller than — e^/(37r^) at second order. 

3. The exponents of the density-density and current-current response functions at zero transferred momentum 
(third and fifth channel in the left table) are vanishing. Actually, using the second WI in Eq.(3.13) it can 
be easily proved that these two exponents are zero at all orders in renormalized perturbation theory. This is 
indeed a necessary prerequisite for having finite conductivity and a semi-metallic behavior also in the presence 
of interactions, as expected [49]. 

8. RENORMALIZATION GROUP ANALYSIS IN PRESENCE OF A MASS TERM 



As discussed in the introduction, it is natural to investigate the effects of a mass term coupled to the local order 
parameters associated to the largest critical exponents. In this section we discuss how the RG construction is modified 
by the presence of a mass. Consider, e.g., the Hamiltonian in Eq.(2.27) (similar considerations are valid for CDW, 
AF masses or for the "Haldane mass", i.e., the mass associated to the local order parameter in the fifth line of the 
right part of Table I). The model is still invariant under the same symmetries (l)-(8) described in Section 4B, the 
only novelty being that now the value of jo appearing in the very definition of the Kekule mass, see Eq.(2.27), should 
also be changed under the symmetry transformations; i.e., 

• under (4), jo ^ jo - 1; 

• under (6.b), jo -^ r^jo; 

• jo is left invariant in all other cases. 

In the presence of the mass term, the multiscale integration of the generating functional is performed in a way very 
similar to the one described in Section 4 for the massless case. After the integration of the ultraviolet degrees of freedom 
and the definition of the quasi-particle fields, one immediately realizes that the presence of the mass produces new 
local relevant terms of the form ^^ ^ ^w"/! "k~o- w ^k~cr -w ~^ ^■'^■^ which can be step by step inserted into the definition 
of the fermionic propagator. The symmetries of the model imply the following conditions on the complex constants 
^(io). 



uj^h ' 



1. using (4), we find that A 



Oo) 
uj.h 



^A 



Oo-i). 
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2. using (5), we find tliat A^^^^^ = [^^lH.hT 

3. using (6.b), we find tliat A^^"^' = A^''^^^^ 



This fixes tlie structure of tfie effective mass term, in tire form: 



£W vWfxi, A) = Af — V '|(^'')+r- . i^iv''^^ (8 1) 



k' 



cj.cr 



witli A^ a real constant suclr tliat A^ — Aq is tlie same as tfie one in Eq.(2.27). Tire term in Eq.(8.1) is the only 
extra contribution to the local part of the effective action. In fact, the localization procedure can be now modified 
by requiring that the localization operator extracts from the kernels with 2n + m + p = 2 (or = 3) the first two (or 
the first) terms of a Taylor series in (k, p,q) and in {A^}fe<o: this is because all the terms in perturbations theory 
proportional to the mass itself have a dimensional gain of (9(2(i-'=°">'*-'= )'») with respect to their "massless" bound, see 
[17, Sections 3.2 and 3.3] or [48, Section 5.3.2]. Therefore, with this new definition of localization, the only possible 
extra local term in the effective action is a quadratic bilinear in the fermionic fields as the one in Eq.(8.1). Note also 
that such term can be rewritten in rclativistic form as: 

>C(faLv('')(*,A) = -Af-i,^V^g^;+[e^^.o7.^3]V,g^)- , (8.2) 

with0,„ = f (jo-1). 

Remark. The masses corresponding to the CDW, AF and Haldane mass terms can be worked out along the same 
lines. It turns out that these effective masses have the form: 

C^L^sY'V^'K^. A) = A$r^>^^ Y. *^,S>3*gt!; , (8-3) 

k',CJ,cr 

C'r^Llv'^-^^.A) . A^^ E -^g2>3*gt!; , (8-4) 

k'.W.cr 

c^^Lv^'\^,A) = A-^ E -*g':!:-3*g2; , (8.5) 

k',a;,cr 

which can be rewritten in rclativistic form as 

>cL^aT)v('^)(*,A) = .A^^^-L^^^^^;+^3^(^^)- , (8.6) 

' k',CT 

C\^LlV^^\^, A) = ^Ar -^ E -^S'^3^g^ , (8.7) 

P k',<T 

/:L^aLv('^)(*,A) .= ^Af-i^E^gt'^375^gr , (8.8) 

' k',<j 

which are the same as those considered in [26, 30]. 

As mentioned above, the fermionic bilinear in Eq.(8.1) is inserted step by step into the definition of the fermionic 
propagator. As a consequence, the effective propagator at scale h acquires a mass gap of size A^. Of course, this 
gap is not visible as long as jA^j ^ 2^ . Therefore, if the bare mass Aq is small, the multiscale integration and the 
dimensional bounds remain unchanged up to a scale ft-o such that A^ ~ 2''". At that point the propagator has a 
mass of size comparable with 2'*° itself, and we can integrate the fermionic degrees of freedom associated to scales 
< /lo in a single step. From that on, we are left with a purely bosonic theory. The symmetries of the theory can be 
still used to prove that the only local terms in the effective action for such a bosonic theory are photon mass terms 
of the form 'Yliu'^^^t^-hAiT A\r , h<0, with D^h independent of jg, see [48, Section 5.3.2]. Since the Kekule mass 
term in the Hamiltonian does not break gauge invariance, the same argument used in Section 5 A to control the flow 
of i^fi^h can be repeated here to show that i>^ji — O(e^), which allows one to safely integrate all scales up to — oo. 
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We are left with the problem of computing the scale ho that separates the massless and massive regimes. In order 
to do this, we need to control the flow of A^, which is driven, as usual, by a beta function equation: 



Ati_ Z, 



'i-l _ '* /'l I flA 



(l+/3^.). (8.9) 



Proceeding once again as in Appendix D, we find that 



The resulting flow is: 



I^K,h = ?2 log2 + 0(e4) + 0{e^2^) + 0(6^(1 - y)2^^onst.)e^h-^ (g^^Q) 



Af^AoM2-''-\ with rj^^r^^-r^^^ + Oie^) (8.II) 



and Ao(w) = Ao(l + O(e^) + 0(1 — v)). Therefore, the equation for the dressed electron mass A = 2^" becomes: 

A = 2'"' = Ao(w)2-'''^''« ^ A = [Ao(w)] '/('+''^^ , (8.12) 

which proves Eq.(2.28). 

9. GAP EQUATION 

The variational equation corresponding to the minimization problem Eq.(2.29) is: 

f 

K 



^. -'-{&, (9.1) 



where (•) is the statistical average in the presence of the phonon field (ps,j. We now want to check that the distortion 
(j>g- defined in Eq.(2.30) is a stationary point of the total energy, provided 0o ^nd Aq are chosen properly. The 

stationarity condition Eq.(9.1) with (j)gj ~ 4'sj in the limit /3,i — > oo is equivalent to the two following coupled 
self-consistent equations for (f>o and Aq: 

'^o-^Z^ nl i (27r)3 27r|Sr^^-^ ^ ^■)\p ^i.+p,.\a) e + c.c. , [M.Z) 

n>0 J \ J II 

Aq _ g^ >p (^e)" r dp' dk' r , I .1 H*" a+ 6" \Uo)„-^k'iS,-S^)-^p-{S,„-S^). 

3 - « Z^ n\ J {2nf27T\B\^^^-^' ^'^Jp' V+p'+p-.A'+P^^,'^/ ^ ^ + ^-^^ ' 

n>0 ./ V ; II 

where [rj^iSj ■ A.)] was defined in Eq.(B.2), (^■)^-'"' is the statistical average in the presence of the phonon field 

't^Si ^^^ ^^ denoted by a^^,b^^ the first and second components of the spinor ^^o-' respectively. These equations 
are well defined provided that the right hand sides of the two equations are independent of j, jo and to. Using the 
symmetries of the model, we find that (see Appendix C 3 for a proof) 

(i<Xj^"' :=E^ f7^i[vHSr^-)r:-UpAJ'"'=mAik) , (9.3) 



n>0 



a.' 



k +Pjj,o-''J A^ ^[ J (27r)3 L '• ^ J 'Jp' k'+p'+py'^,o- k +pp,<T/ \ J \ J 



ri>0 -^ ^ ' 

with v4.(k) and S(k') two functions, independent of cj and jo, transforming as follows under the discrete symmetries 
of the model: 

A{k) = A{Tk) = A{Ik) = A*{Rhk) = A{R^k) , (9.4) 

B(k') = B(Tk') = B(/k') = B*{Rhk') = B{R^k!) , (9.5) 
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The dimensional bounds following from the multiscale integration, combined with an explicit lowest order computation, 
show that, if |k'| ~ 2'* with h> Hq (here /ig is the scale defined in Eq.(8.12)), 



Zh kl + vl\n{k' + j^p)\^^ 

3 ^h A;2 + i>2|f](fc'+p^)|s 



^^'.,,,X,,,jf' = ^ ,. .Zr.'L.A ^ + ^'(k')) , (9.6) 



where the correction terms j4'(k') and i?'(k') satisfying the same symmetry properties as j4(k) and -B(k') in Eqs.(9.4)- 
(9.5) and of order C>(e2) + 0(2'') + 0(A,,2-''). If we plug Eqs.(9.6)-(9.7) into Eq.(9.2), we are led to the self-consistent 
equations 



" « t"^A<|k'|<i2^l'8|^(k')fc2+^;(k')2|r!(fc'+p^)|2 ' 

o2 /• dk' 1 A(k') , , 

■ ' ^ ' (9.9) 



« yA<|k'|<i27r|S|Z(k')fc2+^,(k')2|^(fc.+p^)|2 

where the "~" indicates that we are neglecting higher order corrections and, for k' small, Z{\<l') ^ |k'|~'', A(k') ^ 
Ao|k'|~''^ and 1 — w(k') ~ (1 — w)|k'|''. It is apparent that Eq.(9.8) can be solved by fixing 0o to be a suitable 
(positive) constant of order g^ . 

On the other hand, Eq.(9.9) is equivalent to Eq.(2.31) and leads to Eq.(2.32) and to the conclusions spelled out 
after Eq.(2.32). In particular, at small coupling, the critical phonon coupling gc above which the gap equation 
admits a non trivial solution scales like gc ~ v^- This can be proved by observing that the range of p such that 
the denominator of the integrand in the r.h.s. of Eq.(2.32) is larger than (say) |w is contained in the interval [0, p*], 
where 1 — (1 — u)(p*)'' = |w, which gives, for small v, an exponentially small p* , i.e., p* ~ e~'"'^'^^'] in other words, 
at small coupling, the range of momenta such that the effective Fermi velocity is substantially different from the bare 
one is exponentially small, and gives no relevant contribution to the r.h.s. of Eq.(2.32). However, the gap equation 
can be naturally extrapolated at intermediate to strong coupling. As noted in the introduction, the larger the charge, 
the smaller gc\ possibly, gc goes to zero for large enough values of the electric charge. 

Acknowledgements. A.G. and V.M. acknowledge financial support from the ERG Starting Grant CoMBoS- 
239694. 



Appendix A: Functional integral representation and gauge invariance 

In this Appendix we sketch the derivation of the functional integral representation presented in Section 3 and prove 
its gauge invariance properties. 



1. Functional integral representation 

Using standard methods of many body theory [44] and the explicit form of the electron and photon propagators 
discussed in Section 2 A, we can write the partition function Tre^*^^^ of our model in the Coulomb gauge as 

Trle-/^^-} i' Pi^d^)p(C)^dA)e^i.^'^\ (A.l) 



Tr{e-'3(^X+-f^A)} 

which should be understood as an identity between power series in the electric charge e. In Eq.(A.l), P(d^) is the 
same Grassmann gaussian integration as in Eq.(3.4), P^'~^\dA) is the gaussian measure 

p(C){dA) = ^^ n \{d^cM.„dlTnM.„c^^{-,^ Y. 4p[*^'''(p)].^-'4-p}, (A.2) 

peV+ , «=1.2 pGX'fl.t 
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with wl, (p) the covariance matrix in Eq.(2.14) and, if and ipppi{x) ~ 5{xo)ip{x + [p — p')5i), 

cr,(y'=-\i 

The quartic fermionic term in V can be ehminated by a Hubbard-Stratonovich transformation, at the cost of intro- 
ducing an extra component Aq^x of the photon field. This yields a representation of the partition function in terms 
of a functional integral analogous to Eq.(3.1): 

P(d*)p(^)(di')e^(*^'^) = /"p(d7A)P«(dA)e^(*^^)| , (A.4) 

where V(^,j4) is given by Eq.(3.6) and P^{dA) by Eq.(3.5), with w^'*^ replaced by w^ := 'w^'^°° . The functional 
integral representation for the observables can be derived along the same lines. 

2. Derivation of The Ward Identities 

The derivation of the Ward Identities involves manipulation of the functional integral. Therefore, it is crucial to 
define with care the cutoffs used to regularize the functional integral and how they are eliminated. First of all, keeping 
the IR cutoff at scale 2'* fixed, we introduce an ultraviolet cut-off, by replacing the fermionic and bosonic propagator 
by XM(fco)'S'o(k) and Xm{po)wJ^'^ (p), with XM(fco) = x(2^^^|A:o|)- Next, we replace the compact support cut-off 
functions by new functions, exponentially close to them but with full support in R; i.e., we replace x(t) by a smooth 
function x'^i't) such that x'^{t) = x(0 = 1 for < i < i and x(i) < x'^i^) < x(0 + ^e^* for i > i. Finally, we let the 
Grassmann field live on a finite lattice, both in the space and time variables: pick iV G N and define: 

A^^^={xo = ^: < n < TV} X A , ^"^ •= {ko = j{n+^): < n < N} x Bl , 

Let ^ be a Grassmann field on A^ with antiperiodic boundary conditions in the xo variable, with propagator 

Soika,k) :— S'o(dAr(fco), fc), where (iAr(fco) — i- — jTff—^- We shall think of Eq.(3.1) as the M,N ^ oo and e ^ 
limit of a regularized functional, with the limits taken in the following order: 

gW«^''-($,j,A)^ lim lim hm e'^'^'^^'^^}, ^w^^%(^,J.x) ^ / PM,iv(rf^)P^'* (dA)e'^(*'^+^)+^(*'^+'^'*)+(^^*) , 

M— !-oo e—>^) N—>-oo J 

(A.5) 
where the interaction and source terms are the same as in Section 3, with the only differences that J dx should be 
understood as a shorthand for -^ ^ ca*"' ^^^ ^^'^ term in the second line of Eq.(3.6) should be replaced by 






'/3 

^+ ^^- .. U^^ii'" Af,,i^^^^,,s) _\\j^^+ qy- /gie/o ^o,(xo + = ,x+ii) _ 1^ 



Now, the regularized functional e '^'^."^ ' ' ' is a non-singular finite dimensional integral. Therefore, we can freely 
perform unitary change of variables of the fields ^, without affecting the value of the integral. In particular, by 
performing the unitary phase transformation ^^ „ „ — > ^^ ^ e^'^'^°^y>^+{p^'^)°i-) ^ we find: 

0=T^Wi^"^(<i>,J + aa,Ae*^") -f A|,_jv(P; *, ^, A) , (A.6) 



aAp'^'^' 



where 



Ai.,^(p;$, J,A) = J- log /pM,^(d*)Plj''*(ciA)e'^(*^^+'^)+«(*^-^+^^*)+(^^*)+'^(*^*) , (A.7) 

0(Pp J <t>=o 
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and 



C(*:<^):-«27T^ E ^*k+p..a/,iv(k,p)*k,.. 



P^L^\Sl\ 



(JV) 



kee^.z. 



,-1 



, -1 



Cl^,jv(k,p)= [(X|^(fc0+P0)) -lJSAr(k + p)ro(p)-[(x|^(fco)) -lJro(p)i?Ar(k), (A.8) 

where -BAr(k) = [S'o((iiv(fco), k)\~^ . The correction term A^ ^ is due to the presence of the imaginary time UV cut-off 

in PM,N{d^), which shghtly breaks gauge invariance. It can be studied by a muffiscale analyisis of the functional 
integral in Eq.(A.7), using ideas similar to (but much simpler than) those of [18, Appendix D] (using methods first 
developed in [4, 5]), where the correction to the Wis due to a fixed fermionic UV cutoff at momenta of order 1 were 
studied. Here the main difference is that the UV cutoff is only on the imaginary time coordinates and it cuts off 
momenta of scale 2*^ . By the support properties of the kernel C%j ^ (k, p) all the contributions to A^ ^ (k, p) have at 
least one loop momentum flowing at scale 2^^ . Moreover, since the scaling dimensions in the UV theory for the time 
variables are all negative, see [15, Appendix C] and [48], it is easy to show that all such contributions go exponentially 
to zero as M, N -^ oo: 

lim lim lim A^^ ^(p; $, J, A) = , (A.9) 

for all fixed p, which proves Eq.(3.11). More details on the proof of Eq.(A.9) can be found in [48]. 

3. Independence of the functional integral on the choice of ^ 

In this subsection we show that the averages of gauge covariant operators do not depend on the choice of the gauge 
fixing parameter ^ appearing in the bosonic integration measure. This allows us to work in the technically convenient 
Feynman gauge, corresponding to the choice ^ = 0. 

Let F('^, A) be a gauge invariant function, i.e., F('^, A) = F{'^e^'^°', A + da). We want to show that 

JPid^)pi'^'{dA)e^(^^^^)Fi^,A) _ 
^ JP{d-^)Pi-h'{dA)eVi-^^^) ^ ' 

As in the derivation of the Wis, the proof requires a number of manipulation of the functional integral. In principle, 
we should proceed as in the previous subsection, by first introducing proper UV cutoffs in the time variable, then 
keep track of the possible correction terms produced by the manipulations and finally discuss the vanishing of the 
corrections in the limit where the UV cutoffs are removed. However, the result is that, once again, these correction 
terms vanish in the limit, so here we will neglect them from the very first moment. 
With obvious notation, let us denote the gaussian measure in Eq.(3.5) by 

2Me-5('4,(-«'''*)-^A) 
P«''' (dA) = ^^^^^-^ , (A.ll) 

and, given a function of ^, A, let us define 

E^, (0(*, A)) =. l^^^ ,,^, ..Ml. • (A.12) 

Note that, letting Oi := e^ F and O2 '■= e^ , Eq.(A.lO) can be rewritten as 

rp(dvi/)Ei (©!(*, A)) 



■JP{d^)Ei{0,{^,A)) 



= , (A.13) 



where Oi and O2 are both gauge invariant functions. So in order to prove Eq.(A.lO) it is enough to prove that 
/P((i\l/)5jE^(0(^, A)) = 0, with O a gauge invariant function. To this purpose, note that 

d^Ei(0{^,A)) = -^^ ^ d^[wi^'^'{p)]-l[Ei{A,^pA,.,^pO{^,A))-wlf{p)Ei{0{^,A))]. (A.14) 

M, 1^=0, 1,2 
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Integrating by parts the first term in square brackets, we find: 

1 



d^Ei{0{^,A)) 



2I3\Sl 



J2 ^^a (p)9« [*«^'^* (P)];>«;.f (p)iEl ( 



-i,;,«,'«V.w«7^!^(!^^ 



fi^UjUi ,1/2—0,1,2 



c* j4jy J ^ p a A^2 , — p 



/J, 1^=0, 1,2 

Using the exphcit form of wj^'^ (p), Eq.(3.3), the last expression can be further rewritten as 



(A.15) 



d^E^JOi^,A)) 



1 V MP3 x(bl)-x(2-'' IpI) 

M, 1^=0, 1,2 



(A.16) 



• [PuPi' - PoiPfj^n,, + p^Ufj^)] 



&" 
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E5i(0(vl/,A + J)) 



,7=0 



On the other hand, by the gauge invariance of O and performing the unitary phase transformation 'i^ ^^ „ 
^x o■p6**'^"'''^'■''^^''^^^ ^^ i^ the previous subsection, we find 



P(d*)E«i(0(*,A)) = / P(d*)E5i(C'(*,A + 9a)) 



which implies 
d 



9dr 






= 



,7=0 



(A.17) 



p=0,l,2' 

Integrating Eq.(A.16) with respect to j P{<M>) and using Eq.(A.17) finally gives the desired cancellation. 



Appendix B: Symmetry transformations 

In this Appendix we prove that both the gaussian integrations and the effective potentials are invariant under the 
symmetries (l)-(8) listed in Section 4B. As already done in Section 4, we restrict for simplicity to the case that 
J = A = and that all the external fields $" but the one with a = K are set to 0. The invariance of the effective 
potentials on scale h and of the single-scale integrations under the stated transformations follows from the fact that 
the bare interaction and source terms, Eqs.(3.6)-(3.7), and all the single scale integrations P((i^(''^)P((iA('')) are 
separately invariant under the same transformations. Moreover, since the single scale integrations are obtained from 
the bare integration P((i^)P''''' (dA) by recursively including the local terms ^C^V''*-' (which is invariant under the 
symmetries (l)-(8), see Eq.(4.21) and [18, Proof of Lemma 1]), one can immediately convince oneself that the desired 
invariance properties of the effective potentials at all scales follow from the invariance of the bare interaction, bare 
source term and bare integration under the analogue of the symmetries (l)-(8) written in terms of the fields ^j. ^ 

(rather then in terms of the quasi-particle fields ^^, ^ ^, as done in Section 4 B). The symmetries (l)-(8), if rewritten 

in terms of ^^cr' '"^^'i ^^ follows. 

(1) Spin flip: ^^ ^ -^ ^^ _^, and A^p, $^ are left invariant; 

(2) Global U{1): *|.^ -^ e*^"-^^_^, with a^- G M independent of k, and A^^p, ^j^p are left invariant; 



(3) Spin S0(2) : rM 
V*k,i,, 



e 



with 9 independent of k, and A^^p, 'i'f„ are left invariant 
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(4) Discrete spatial rotations: if Tk = (fco, e~'^^'^^k) and n_ = (1 — <T3)/2, 

*k,o- ^ "^ *Tk,o- ' *k,o- ^ ^Tk,<T^ ' ^p ^ -L ^Tp , Vj^p -^ 'i'j+l,Tp 1 

(5) Complex conjugation: if c is a generic constant appearing in the effective potentials or in the gaussian integra- 
tions: 

c^ c* ^f ^ ^"^ , A„ ^ -A „ 4^ -^ <i>^ 

(6. a) Horizontal reflections : if i?ft,k = (fco, — fci, ^2) and rhl = 1, rh2 = 3, rh3 = 2, 
(6.b) Vertical reflections : if _R„k = (fco, fci, — fc2) and r„f = f , r„2 = 3, r„3 = 2, 

(7) Particle-hole : if Pk = (fco, -fc), 

*L. ^ **pL. . ^P ^ P^-Pp , *fp ^ *jy-Pp , 

(8) Time-reversal : if /k = (—fco, fc), 

*k,. ^ -*^3*7k,a . *t. ^ -**lk,a'^3 , ip^/i/p, dj^p ^ <ljy,p . 

Now, the proof of the fact that P{d'^) is invariant under (l)-(8) has already been discussed in [15, Section 3.1]. The 
invariancc of P'^'^ (dA) under (l)-(8) is obvious, and so is the invariance of V(^, A) and ,8(^,^4,$) under (1) (3). 
Therefore, we are left with proving the invariance of V{'^,A) and B{'^,A,^) under (4)-(8). As a preliminary step, 
let us rewrite the interaction and source term in momentum space: 



' ' ' k,p,o- j = l n>l \ V 



g-Jk((5j-5i) 



-ieA, 



_;^)ngj(k+p)(5,-5i) Q / k,<T 

o,p*,r+p,.ro(p)*k,.} , (B.i) 

^ ' ' k,p,q,CTj = l n>0 ■ \V -^y t^ " / 

where, if ry^ = (1 — e^*P^j)/(ipJj), 

(Pl'^il) pi+...+p„=p 

and we recall that ro(p) = _ ,, .In the second of Eq.(B.l), [rj"! (Sj ■ A.)] should be interpreted as equal to 

/3|5i|(5p 0- Let us neglect the spin index, which plays no role in the following, and let us denote by 



(*) ^= ^ E E E ^ [^' (^^- • ^■)] 7^i+p . iy.,.(k+p)(.,-..) *k > 

k,pi=ln>l ' \\ ^J '^ I 



(B.3) 



the term in the first line of Eq.(B.l) and, similarly, 

(**) := -le J2 io,P*k+pro(p)*k , (B.4) 



k,p 

3 



iK)---EE^f.E^-^b^i^.-^X^H,pJ,_,,' '""t"^>k-- (B.5) 



n\ - — -I' ■- ' t- ■ -1 y / -[^\ngi(k+p)(5j-5i) Q 



k,p,qj — 1 n>0 

Let us prove that these terms are separately invariant under the symmetries (4)-(8). 
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Symmetry (4-)- The term (*) is changed under (4) as 



(je)" - --^ .^_.„., / n p-ik((53-(53) 



k,p j — 1 n>l \ ^ '^ / 

where we used that 77^ ~ rjijt . Using the fact that [77:*^ ((5^ -(T^-^yl.))] ™ = [77:'^ ((5j+i-j4.)]™ andp-^j — (Tp)-5j+i, 
we see that the r.h.s. of Eq.(B.6) can be rewritten as 

k.p j = l n>l \ ^ -^ / 

which is the same as (*), as apparent by the change of variables (Tk, Tp) -^ (k, p) in the sum. The proof of the 
invariance of the term {K) under (4) is exactly the same. Regarding the term (**), it is changed as 

(**) -^ -*e E ^o.^P*T(k+p) ( ^-Il„3 ) *Tk,. . (B.8) 

k,p \ ' / 

which is the same as (**), simply because /> ~ ro(Tp) and we can perform the change of variables 

(Tk, Tp) -^ (k, p) in the sum. 

Symmetry (5). The term (*) is changed as 

w-^EEE^f [^^(^vM)r;*VpL_, °.pK..-., , pik, (B.9) 



k,p j — 1 n>l 



(—\\np-i{^+p){Sj-5i) 



where we used that (?;i)* = tj-' . The r.h.s. of Eq.(B.9) is the same as (*), as apparent by the change of variables 
(k, p) -^ (— k, — p). The proof of the invariance of the term [K) under (5) is exactly the same. Moreover, the term 
(**) is changed as 

(**) ^ +ze^(-io,-p)*+k-pro(p)*:k > (B-10) 

k,p 

which is the same as (**), because rQ(p) = ro(— p) and we can perform the change of variables (k, p) -^ (— k, — p) in 
the sum. 

Symmetry (6. a). Using the fact that RhSj = —Sr,j and 77^ = e^^^^^rj^j^-' , the term (*) is changed as 

k.pj = ln>l ' \ V ^^ "= " / 

(B.ll) 
which can be rewritten as 

^EEE^[^^^^('5'--^^)];:p*^.(k+p)f, ,.„^_.(k^^^ '"'7'^')*i..k, (B.12) 

k.p j = l n>l \^ ' / 

which is the same as (*), simply because e'^^^^^-^i) = g-^RhH^r^j-Si) ^^^ ^^ ^^^ perform the change of variables 
(i?hk, RfiP) -^ (k, p) in the sum. The proof of the invariance of the term (K) under (5) is exactly the same. Moreover, 
the term (**) is changed as 

(**) ^ -ze^io,i?.pe*P'^*t(k+p)'^iro(p)'^i*]^,k : (B-13) 

k,p 

which is the same as (**), because CTiro(p)cri — e~*P^iro(i?/ip) and we can perform the change of variables 
(i?ftk, Rhp) ^- (k, p) in the sum. 
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Symmetry (6.h). Using the fact that R^8j = 8r^j and 77^ — ?7^""'p, the term (*) is changed as 

k,p j=l n>l ■ \ V J^7 "' " / 

which is the same as (*), simply because ^-^'^'y^i-^^) — ^-iR-c\i.{Sr^i-&i) ^-^^^ -^g ^g^j^ perform the change of variables 
(i?t,k, RvV) ^ (k, p) in the sum. The proof of the invariance of the term (K) under (5) is exactly the same. Moreover, 
the term (**) is changed as 

-le J2 ^.i^.P*l(k+p)ro(p)*«„k . (B.15) 

k,p 

which is the same as (**), because ro(p) = ToiRyp) and we can perform the change of variables (i?^k, _Rt,p) -^ (k,p) 
in the sum. 

Sym,metry (7). The term (*) changes as 

•^ fie)" / g-ik(6j-5i)\ 

(*)^*EEE^[^'M^-^)]-pp*pkL . ,(k+p)(,,-,,) *p(k+p) = 

k,pi=ln>l ■ \\ ^J '^ " / 

= *EEE^[^^('5.-^)]-Pp*pk , ,,„,™,-.,) n r^C'+p)' (^-^^^ 



k,p J — 1 n>l 



(— l)"e*^'^^^^~^i' 



which is equal to (*), as apparent after the change of variables (P(k + p), — Pp) — )■ (k, p) in the sum. The proof of 
the invariance of the term (K) under (5) is exactly the same. Moreover, the term (**) is changed as 

(**) ^ -te ^ io,-Pp*+kr;^(p)*p(k+p) ' (B-17) 

k,p 

which is the same as (**), because r^(p) = Fofp) = Fof— Pp) and we can perform the change of variables (P(k + 
p), — Pp) -^ (k, p) in the sum. 

Symmetry (8). The term (*) changes as 

(*) ^ -*EEE ^[^'('5^- • ^oiIp^W)-^^ i)ng.(k+p)(.,-6,) j"'*^^'- = 



k.p j—1 n>l 

^EEE^[^^(^.-^^)];:*W)i. iyv('^"+p)(..-M ' V '')*-' ^^-^'^ 

k.p 7 = 1 n>l \ ^ ' / 



""(*e)",,. ^-.,.*„^, / e-k(^.-^i)^ 



which is equal to (*), because e^^'^^^^-'^i) — g-«/k(5j-(5i) ^^^ .^^ ^^^-^ perform the change of variables (/k, /p) — > (k, p) 
in the sum. The proof of the invariance of the term [K] under (5) is exactly the same. Moreover, the term (**) is 
changed as 

(**) ^ -le Y, ^,/p'I'|(k+p)'^3^o(p)a3*7k ' (B-19) 

k,p 

which is the same as (**), because <J^TQ{\i)(j^ — Tofp) = Fof/p) and we can perform the change of variables 
(Ik, /p) -^ (k, p) in the sum. This concludes the proof of the invariance properties stated in Section 4B. 

Appendix C: Symmetry properties of the kernels 

In this Appendix we exploit the lattice symmetries (l)-(8) listed in Section 4B to prove the invariance properties 
Eqs.(4.20)-(4.21) of the local terms in the effective action stated in Section 4C. We will start with the "relevant" 
terms with 2n + m + p = 2 (i.e., the terms of the form AA, ^^ A or ^+^1'^) and we will then proceed with the 
"marginal" terms with 2n + m + p= S (i.e., the terms of the form A'^^'^', <I>-^v]/+\]/-^ AAA, ^^ AA or ^^^^A). In 
the following, we shall drop all the unnecessary labels (including the scale and spin labels), to avoid an overwhelming 
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notation. We will think the operators T, Ry, Rh, P, I appearing in the symmetry transformations as 3 x 3 matrices 
acting on the /x-indices, with 

/l 

T = cos ^ - sin ^ I , Rh = \ -I \ , Rv ^ \ 1 \ , P = -/ = | -1 | . (C.l) 

\0 sinf^ cos^ 

1. The "relevant" terms 

The structure of the local terms of the form ^+^1'^ has already been studied in [15, Lemma 3], where the first 
equation of Eq.(4.21) was proved. Let us then look at the terms of the form AA or $^A. 

The AA kernels. The contribution to the effective potential quadratic in the A field is proportional to 

^ipVK(p)i_p, (C.2) 

p 

for a suitable 3x3 matrix-valued kernel W^(p). Imposing the invariance of Eq.(C.2) under symmetries (4)-(8), we 
find: 

Wip) ^^ TWiT-^p)T-^ ^^ [Wi~p)]* ' = ' RhWiRhp)Rh ^ = ^ i?„W^(i?.p)i?. ^^ IW{Ip)I . (C.3) 

Let Wp^(O) =: f^^ and daWfj_,y{0) ~: (^^^q,; the properties (6. a), (6.b) and (8) in Eq.(C.3) imply that i^^^ — 5^j_i^v^, 
v'a " 0, while (5) in Eq.(C.3) gives i/^ = v*' finally, the property (4) in Eq.(C.3) implies that vi — V2. This proves 
the second equation in Eq.(4.21). 

The '^^ A kernels. The contribution to the effective potential of the form $^A is proportional to 

J p 

for a suitable vector-valued kernel W^{p) ~ (M^j^(p), T4^j^(p), M^^^(p)). Using the symmetries (7) and (8) we get: 

Wf{p) ^^ Wf{-Pp)P ^^ Wf{Ip)I , (C.5) 

which implies that W^{p) — 0, because P = —I. 

2. The "marginal" terms 

In this subsection we study the structure of the terms of the form Avp+vj/-^ ^^vp+v]/-^ AAA, ^^ AA or $^$^A. 

The A^+^~ kernels. Since A has an UV cutoff that suppresses modes Ap with p close to ±(pj — p^) (and to 
its images over A*), the only non zero terms of the form A^t'+^I'^ are those with the two fermi fields associated to 
the same w index. [Of course, if we were interested in studying the theory without the UV cutoff on the photon field, 
then terms of the form A^+\I'Z(^ would be allowed, and their structure could be investigated by the same methods 
used here.] The contribution to the effective potential of the form A^^^J is proportional to 

k',p,a; 

for a suitable tensor-valued kernel W[j(k',p) — (IVo_(j(k',p), W'i^(^(k',p), W'2,(^(k',p) (each component M^^^i^(k',p) is 
a 2 X 2 matrix, acting on the p indices of the grassmann fields). The invariance under the symmetries (4)-(8) implies: 

VK„(k',p) '=' e-*(P^+'''+P)(^i-^^)"-[VK^(T~ik',T-ip)T-i]e'(P^+*'')(^i-^^)"- ^=' -#_!^(-k', -p) ^="' (C.7) 

(6-) e-'P''a,[WURhk',RhP)Rh]<Ji '-^ W^-c.(i?.k', i?„p)i?, ^^ W^jP{k' + p),-Pp)P ^^ -a3[WUlk\Ip)I]a3 . 
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Let us define 

3 

^^,..(0, 0) = ^ al^^a, , a:: - K,^, a^,^, aU , (C.8) 
with CTo = 1. The properties (C.7) imply that (summation over repeated indices is imphed): 

a>, 'i' -[a"_JVo-[aL^]Vi + [a^JV2-[a3_JV3, (CIO) 

a^^^cr^ = [a°i?ft](To + [alRh](7i - [alRh]<J2 - [alRh]<J3 , (Cll) 

a>, * = ' K^i?„]a, , (C.12) 

a>, '=' [aO ^P]ao + [aL^P]ai - [a2_^P]a2 + [a^l^Plag . (C.13) 

a>, 'I:' -[a;i/]ao + [ai/]ai + [a2/]a2-[a3/]a3. (C.14) 



We recall that T,Rh,Rv,P and / are the matrices in Eq.(C.l); the notation, e.g., [aT^^] indicates the row 3-vector 
obtained by matrix multiplication of the row 3-vector a times the 3x3 matrix T^^. Properties Eqs.(C.ll) and (C.14) 
imply that 

aJi = (a°,^, 0, 0) , ai = (0, 0, fli,^) , a^ = (0, a^ ,^, 0) , a^ = (0, 0, 0) , (C.15) 

while from Eq.(C.12) we get that 

a^c^ = «o,-^ =: -*^o , al^ = -o-i-^ =: -wAa , a^^^ = al_^ =: -Ai , (C.16) 

with A^ e R, thanks to Eq.(C.lO). Therefore, 

WoAO, 0) = -iXo<Jo = jAoEO , #1,^(0, 0) = -Aia2 = iXiTl , m^O, 0) = -Aaw^i = lAaE^ . (C.17) 

Moreover, thanks to Eq.(C.9): 

e-i^ j + [^^^"] (^ei- o' j' ^^-^'^ 

which gives 02 ^ = wa^ ^, i.e., Ai = A2. This proves the equation in the first line of Eq.(4.20). 

The <i)^v[/+v]>- kernels. The contribution to the effective potential of the form ^^'i^'^^, is proportional to 

E '^,':p'+p^,-p-*k'+p',.^^.,.o(k',P')*k',.' , (C.19) 

k',p 

for a suitable 2x2 matrix-valued potential VF^(k',p'). The symmetry properties (4)-(8) imply that: 

(6.a) 



Wf^(k',p') ^=^ e-^(P-+'^'+P')(^^-^^)"-T4^^i,^(r-ik',r-ip')e*(P-+'^')(^i"^^)"- *=^ [Tyf_^(-k', -p')]* ^ = ' 
(•L") e»(p-.-p-;+p')(^.-^.)^^^^^^.^(i?,k',i?,p')^i ^ = ^ W:'^^,_JRX,R.P') ^^ 
= [W^^-.',-.)(^(k' + P'), -^P')] "^ = -'^3W,^.(/k', /p')a3 , (C.20) 



Let us define W-^^{0, 0) == J2l=a ^j.^^^n^- Then, property (8) in the third line of Eq.(C.20) implies that b^^ = b^^ = 0. 
Moreover, using properties (6.b) and (7) in the second and third lines of Eq.(C.20), 
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which implies, in particular, that b^ ^ ~ and 6j / .^ — b±, for two suitable (w-independent) constants b±. Using 

(5) 

property (5) in Eq.(C.20), b^ ^ ~ {b] _i^)* , we see that both these constants are real, b± G K. Property (4) in 



Eq.(C.20) implies that, iitj — (a;,a;') 



so that, finally. 



(4)/ e-'^(&}_i,^-z62_^j 



"- "- *e-¥(6i +z6? ,J 



W^^.,±.)(0,0) = 6±(^^_^(^._^, ^ j, (C.2f) 

which proves the second of Eq.(4.20). 

The AAA, $^AA and (^'^(^^A kernels. The contribution to the effective potential of the form AAA or <^^<^^A 
are proportional to, respectively, 

Y^ ^Pi,Pi^M2,P2^M,-Pi-P2W^(mi,M2,M3)(P1.P2) Or XI ^J^,Pi*jI,P2^M,-pi-P2W^(fij2),M(Pl'P2) ' 

P1.P2 Pl,P2 

for suitable kernels W^(^i,p2,A'3)(Pi' Ps)' ^(fi^a) u(Pi'P2)- Using the invariance under the symmetry (7) + (8) we find 
that W^(^i,^2,^3)(pi,P2) = -VF(a.i,m2,m3)(P1'P2) and W^(^^ j^)^^(pi, P2) = -W^(7^ j,),^(Pi, P2), that is, they are both 
identically zero. On the contrary, the contribution to the effective potential of the form ^^ AA is proportional to 

Z2 '^J!pi"^Ml.P2^A'2 -Pi-P2W^j\^^i,^2)(Pi'P2) ' 

P1.P2 

5,^1.^2 

for a suitable kernel W'^:^^^^^^^(pi, pa). Its local part #'^^^^^^^^(0,0) =: Aj^^^^^^^j satisfies: 

\^ (1) rp.K rp-l W r,^^]* (6:") n 5,A: n (^i.) p xK p 
Aj — 1 Aj_il — [Aj J — JXhA^^jJXh — nvA^^jIty , 

which imply 

/a 0\ fa Q ^ \ /a 

Af= 60 Af = i(6 + 3c) ^(c-5) , Af - i(6 + 3c) ^(6 - c) 

\0 c/ \0 ^(c-6) i(c + 36); VO ^(^"c) i(c + 36), 

for some real constants a, 6, c. This the general symmetry structure of the terms in the first line of Eq.(4.28). 

3. Symmetry structure of the kernels in the presence of the phonon field 

In this Appendix we prove Eqs.(9.3),(9.4),(9.5). Let us define, for all k, k' 7^ p^, 

-7+ h- ))^'"'^ {{at -Jk' ^^^'"^ 

m ' ''"'^^^^- ' n{k' 

Using the symmetries listed in Section 4B and in Appendix B and proceeding as in Appendix C, we find: 

Af )(k) ^^ A^:t'\Tk} ^^ [A^"\-k)Y ^'=^ [Aii"^{R,.k)Y ^'='^ A^"\R.k) ^^ [A^'"\Pk)Y ^^ A^"\lk} , 

5S(k') ^ S^^:i'(rk') ^ [B%{-k)r '^^ [B^-j^^jR,k)r <^' i?t?i,(i?.k) <5 [B%iPk)Y <5 s(;^«)(/k) 

Using the properties (6.a)+(6.b)+(7) in the latter equation, we find: A^'''\k) = A^'"''^"\\<.) and B^J^J(k) = B^_!'^^f{k). 
Combining this with property (4) we get: 

A^^^(k) = 4^'(Tk) = A'i^Tk) = A'^^\T^k) = 4^^(^^k) , 

Af ) (k) = Af ' (k) = Af (Tk) ^ a(2) (Tk) = 4^) (T^k) , Af ) (k) = A« (Tk) = Af (T^k) , 



^Oo)^j^) _ VV°k,^V<x»j ^ B(/«.)(k') := c-»p"^(^^"-^i^ '^'+P^^ k'+p^//j- ^ ^^^^^^ 
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i(io) 



Oo). 



which imphes that A^ (k) = Aj(k). The analogous equations for B imply that B^"-' (k) — i?j(k). Now, using 
properties (4) + (6.b) gives Ai(k) = A2{Tk) = AsiTR^k) = Ai{T'^R^k) = Ai{T'^k), which implies that ^i(k) = 
Ai(Tk); therefore, Aj{k) is independent of j and transforms as in Eq.(9.4). The same argument holds for B. 

Finally, it is straightforward to check that the previous symmetries also imply that, if k = p^, then {{a^ ^b 



\Oo) 



•<V+p--^k'+p-//j 



Oo) 



^k,cr'^k,cr//j 



0, which concludes the proof of Eq.(9.3) for all k, k'. 



Appendix D: Lowest order computations 



1. The beta function for ^jf ^ 



Here we prove the two identities in Eq.(5.19) (the identity in Eq.(5.18) has been proved in [18, Appendix C]). The 
relativistic part of the second order beta function is given by the sum over /ii, /12, h^ of the graph in Fig. 6, with the 
constraint that min{ft,i, /12, /is} = h and that hi — h Cz {0, 1}. In the graph, the external wavy line with transferred 
momentum p'i/ — p'p represents $^ c^i_ ^,2 and the straight internal lines correspond to Dirac propagators on scales 

hi,h2', ii oji = oj and uj2 = ±0;, then the graph represents a contribution to /^^F^ or /?^FJ , where F^ was defined 
in Eq.(4.20). 



v7 - Pf 




FIG. 6. The lowest order contribution to the beta function for Z'^^'^^ . 

If we neglect corrections of order 0{e'^) and 0(e^(l — v)2'^^ ^) for some c > 0, in the computation of the graph we 
can replace e^^^ by e, fh{^) by //i(k') (see Eq.(4.30) for a definition of fh) and Vh by 1. As a result, we find: 



'^'-•^-'% L /(W^'-'^'^'-'^'^'-'^'^R^-S^F^Tl^^' 



"'J,-^ fO 



1^=0 hi>h 

min{hi}—h 



ipo^±^ +ip-T±, 



-n., (D.i) 



modulo corrections O(e^) and 0(e^(l — ^;)2('^°"^*-^'^ '*). It is convenient to pass to spherical coordinates. Note that, 

if ^,,(|p|) :- fhipf + 3fh{p)\fh+i{p) + 3fhip)fh+i{p?, then J^°° dp p-^Fh{p) = log2. Therefore, after having 
integrated the radial coordinates we are left with: 



2 „i 



r2Tr 



•F;;; ( - iF° cos + iVl sin 9cos(p + iVl sin 6 sin (p) F^_^- ( - iV'^^^ cos + iV]^^ sin 9cos(p + iF|^ sin 6 sin 0) F^^ , 
which implies 



/3.^ = 



e^ log 2 I 
127r2 



K J "' E r:: (r° r±,TL + riF± /L + r^r± ^tL) C 



(D.3) 



Now, if we pick the plus sign, F^ = cti cos9j — uj(J2 sin^j, where 6j = ^{j — 1). Therefore, using the definitions of 
F^inEq.(4.11), 



-pO-p+ pO ^ -p+ 



plp+ pi ^ p+ 

w-^ WJ " — t^ J ' 



p2 p+ p2 p+ 



(D.4) 
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Plugging Eq.(D.4) into Eq.(D.3) we find jif^ = r^^ log 2, as desired. 
Similarly, picking the minus sign in Eq.(D.3), FJ — e^*"^^cri, so that 



' U) UJ,J —U) U.J ' LO U.J —U LJ,J ' U U.J —U U ,J 

3e^ 



Plugging Eq.(D.5) into Eq.(D.3) we find /3^ = ||7log2. This completes the proof of Eq.(5.19). 



2. The beta function for the other renormahzation constants 

A similar computation can be performed for the beta function controlling the flow of the other renormahzation 
constants. In particular, the beta functions for Z^jjiyi^, ^cdt^/i ^nd Z\pf^ and Z^p^^ are given by expressions 
analogous to Eq.(D.3) with F^ replaced by 0-3 and e*'^''^""(T3, respectively. Using the fact that 

^u'^i^u, = ^■i 1 r^CTaF^^ == 0-3 , F^CTaF^^ = 0-3 , 
we find that P(JJJ^r|^ — j^ log 2, modulo subdominant corrections. Similarly, 

T^O AuOjU- ^ -pO AuOjU-^ T^l iuQjii- iuO ^11- „ -pl AuOTn+^ T~^2iu9jn-^ -p2 iuO^n^^ 

^u^ '^3i-i^ = e ' CT3 , i^e ' e ' cr^l _^ ^ e ' +0-3, l^e ' a^i_^ = -e ' +0-3, 

implies that PcDWh ^ i27r^ log 2, modulo subdominant corrections. This proves Eqs.(7.3)-(7.5). The exponents in 
Table I are computed analogously and we will not belabor the details here. 



Appendix E: Lowest order check of the Ward Identities 

In this section we check at lowest order in non-renormalized perturbation theory the validity of the Wis that we 
used in Section 5 to prove the infrared stability of the flows of the effective charge and of the photon mass. 



1. Ward identity for the photon mass 

We start by checking the WI Eq.(5.6) for the photon mass. At lowest order, the graphs contributing to the dressed 
photon mass are depicted in Fig. 3. As we are going to show here, the sum of the two graphs computed at zero 
transferred momentum is exactly vanishing, for all choices of fi, v. 

The first and second order interactions (of the form ^4^+^^ and AA^^'^^) involved in the computation of the 
two graphs are obtained by expanding the interaction V(^,A) in Eqs.(3.6)-(B.l) up to second order in the electric 
charge: 

k,p,CT j = l \ J 



Defining 



k,pi,p2i=l ^ 



2 -J^ -* f g-ik(6j-5i)\ 

r(k,P):-3i:^'p'^^|_g.(k+p)(.,-.o j' 

2 -J^ - ^ ^ f g-ik(5j-5i) 

Fj„(k,pi,p2) = 3 2^^^,'?^.('^.)i('5,)™( ,(k+p,+p,)(,,-,,) ^ 
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and recalling that v = |i, we can rewrite the interaction as V(^, A) = V2(^, A) + O(e^), with 

k,p,cr 

- 2 03L21S 12 Y. Yl ^',Pi^'™,P2*k+pi+p...r,„(k,pi,p2)*k,a- (E.l) 

^ ' ^' k,pi.p2 l,m=1.2 

Since the term in the last line involves only the spatial components of the photon field, the second graph in Fig. 3 
gives a non zero contribution only if both fj, and v are different from zero. 

If /i = i^ = 0, only the first graph survives, which in the /3, L — > oo limit, if computed at p = 0, gives (using the 
fact that ro(0) = 1) 

where S'o(k) was defined in Eq.(2.11) and \B\ ~ 87r^/(3\/3) is the area of the first Brillouin zone. Using the fact that 

^o(k)S'o(k) = -zao^o(k) , (E.3) 

we see that Eq.(E.2) is the integral of a total derivative, which is zero. 

Let us now consider the case where fx = v = 2. The sum of the two graphs in the /3, L — > oo limit, if computed at 
p = 0, gives: 

■ / 2^ Tr{5o(k)r2(k, 0)^o(k)r2(k, 0)} + ^ | ^ Tr{5o(k)r22(k, 0, 0)} . (E.4) 



2 
Using the fact that 

^o(k)«r2(k,0)^o(k) = ia2^o(k) , (E.5) 

and integrating by parts, we can rewrite Eq.(E.4) as 



^ / 2^ Tr{5o(k) [ - z92r2(k, 0) + r22(k, 0, 0)] } , (E.i 



which is zero, simply because the matrix in square brackets is identically zero. This proves that the graphs in Fig. 3 
with fi = ly = 2 cancel out exactly. Using the symmetry (4), we find that the diagonal terms with fi — v — 1 cancel 
out, too. The non-diagonal terms can be treated analogously. 

2. Ward identity for the effective charge 

Let us check at lowest order in non-renormalized perturbation theory the WI for the effective charge, Eqs.(5.8)-(5.9). 
This amounts to check the cancellation of the graphs depicted in Fig. 4. In order to compute these graphs we use the 
bare fermionic propagator S'o(k) and the photon propagator w^'^ (p) —: w^-'^ Hp)Sfiu with IR cutoff on scale h* , see 
Eq.(3.3). As we are going to show, the sum of these six graphs is exactly vanishing; therefore, the dressed charge is 
equal to the bare one at lowest order. Remarkably, this cancellation does not depend on the presence of the infrared 
cutoff on the photons; this fact has been exploited in Section 5 to derive a WI for the effective charge on all IR scales. 
We shall only consider the cases /i = and /i = 2; the renormalization of the charge corresponding to /i = 1 is equal 
to the case /i = 2, thanks to the discrete rotational symmetry (4). 

In order to compute the graphs in Fig. 4 we need the interaction V(^, A) up to third order in e: 

V(*,A)=V2(VI',A) + 1-^J^ ^ ^ ApiAn.p.Afe,P3*,^+p,+p,+p3,,r;„fe(k,pi,p2,P3)*k,. + 0(e4), 

^ ' ^' k,pi l,m.k=1.2 

P2,P3 CT=ti 

(E.7) 
where V2(^, A) was defined in Eq.(E.l) and 

2J-. ^ ^ ^ / g-»k(5,-5i)\ 

r,™fc(k,pi,p2,P3) = -^r;^,r;^,r7^3(<5,),((5,)„(,5,)fe( ^^(k^p^+p^+p^)(,^_,^) ^ j. (E.8) 



43 

a. Case /z = 0. In this case only the first and fifth graphs in Fig. 4 are non- vanishing. Their value in the /3, L — > oo 
limit computed at the Fermi points k — pf, and at transferred momentum p = is 



»e^(*e.)2 J -^u;(^''*)(p){r,(p^ + p, -p)^o(p^ + p)ro(p^ + P, 0)^o(p^ + p)r.(p^, p)} 



+eV(ze.)2 /-^u;(^''*)(p)a,„|r.(k + p,-p)^o(k + p)r,(k,p)} , (E.9) 

where we defined eg = e, e/ = ve for I G {1, 2} and ro(k, p) := — ro(p). Now, in the first line we can rewrite: 

SoiPF + p)MPf + P, 0)S'o(p'^ + p) = tdoSoip'^F + p) , 

which partially cancel with the second line. We are only left with the terms in the second line where the derivative 
dk„ acts on the kernels F^; however, these terms are identically zero, simply because the kernels F^ are by definition 
independent of fco- 

b. Case /i = 2. Here the situation is more complicated, because of the presence of the second, third, fourth and 
sixth graph in Fig. 4, and because the dk2 derivative can now act on the kernels F^. However: 

(i) repeating the same argument used in the case /i = and using Eq.(E.5), we see that the sum of the first and 
fifth graphs is equal to 



'■J2('^-)^J-(^^^-'"'\p){9k2^-i^ + P' -P)'^0(k + p)F,(k, p) + F,(k + p, -p)^o(k + p)afe,F,(k, P)}^^^^ 

(E.IO) 
(ii) combining Eq.(E.lO) with the contributions from the second and third graphs, we get 



eE(*^-)'/7^^^-''*^(P){[^''2r,(k + p,-p)+zF2.(k + p,0,-p)]^o(k + p)r„(k,p) 



u=0 



+r,(k + p,-p)^o(k + p)[afe,F,(k,p) + jr2.(k,o,p)]| , (e.u) 

which is zero, simply because the terms in square brackets are identically zero, as one can easily check, 
(iii) the sum of the fourth and sixth graphs in Fig. 4 gives 

-^"'^ f 7^^^-"'Hp)[^k,^ui^,P,-p) + ^rMK0,P,-p)]^^^. , (E.12) 

which is zero, simply because the term in square brackets are identically zero, as one can easily check. 

The case /i = 1 can be obtained from the case /i = 2 by the discrete rotational symmetry (4), so this concludes our 
lowest order check of the Ward identities. 
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